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EXECUTIVE  SUMMARY 


This  Final  Report,  issued  in  two  Volumes,  is  a  self-contained 
documentation  of  the  overall  activity  that  Raytheon  has  performed 
under  the  terms  of  Contract  F49620-87-C-0050.  Volume  I  illustrates 
the  program  accomplishments  in  the  most  important,  conclusive  period 
of  contract  performance,  from  August  1988  through  January  1989. 
Volume  II  is  a  collection  of  five  Quarterly  Reports,  published  by 
Raytheon  to  illustrate  project  activity  between  August  1987  and  August 
1988. 


These  past  six  months  have  been  of  great  importance  for  the  project, 
because  of  the  availability  of  tritium  sources  of  neutrinos  at  LANL, 

Los  Alamos,  NM.  This  has  made  it  possible  to  acquire  data  of 

scientific  value,  and  to  answer  some  of  the  questions  that  had 
motivated  the  conduct  of  the  investigation.  These  data  were  processed 
and  analyzed  at  LANL,  in  real  time,  making  it  possible  to  draw 
immediate  conclusions  about  the  meaning  of  the  measurements.  The 
^observations  were  performed  by  using  Prof.  Joe  Weber's  torsion 

balance,  a  room- temperature  instrument  that  was  constructed  by 

University  of  Maryland  under  a  subcontract  from  Raytheon,  and  was 
installed  at  LANL  in  Summer  1988.  The  torsion  balance  was  mounted  at 
a  fixed  location,  close  to  the  edge  of  a  rotating  table  (1  RPM 
rotational  speed)  that  Raytheon  had  constructed  and  moved  to  Los 
Alamos,  NM.  As  the  table  rotated,  the  tritium-filled  container 
(neutrino  source)  and  the  deuterium-filled  container  (that  provided  a 
/*newtonian  forced  reference)  were  sensed  by  the  instrument.  At  the 
/time  of  writing  of\this  Report,  the  results  of  the  measurements  are 
not  fully  conclusive.  A  six-month  Contract  extension,  expected  to 
last  until  31  December  1989,  will  provide  the  final  answer  whether  or 
not  we  could  observe  repulsion  forces,  attributable  to  neutrino 
pressure,  with  the  torsion  balance. 
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The  experiment  consisted,  basically,  of  a  comparison  between  the 
following  two  functions:  (1)  output  of  the  torsion  balance, 
integrated  for  168  hours,  with  a  deuterium-filled  sphere  mounted  on 
the  table  (deuterium  is  non-radioactive,  and  the  only  forces  applied 
to  the  torsion  balance  were  newtonian  attractions);  and,  (2)  output  of 
the  same  instrument,  again  integrated  for  168  hours,  with  a 
tritium-filled  sphere  now  mounted  on  the  table  (providing  a  newtonian 
attraction  force  because  of  its  mass,  with  this  force  decreased  by  a 
repulsion  force  due  to  the  hypothesized  neutrino  radiation  pressure) . 
The  quality  of  the  measurements  was  affected  negatively  by  the 
following  causes: 

(1)  The  torsion  balance  is  characterized  by  a  very  small  damping 
factor  (about  6%),  while  it  would  be  advisable  to  operate  with 
critical  damping  (a  factor  of  100%).  As  a  consequence,  once  excited, 
the  oscillation  train  lasts  well  beyond  the  60-second  period  of  table 
rotation,  and  interferes  with  the  signal  of  the  next  table  rotation. 

(2)  The  torsion  balance  has  an  oscillation  period  (about  14 
seconds)  dependent  upon  room  temperature.  We  noticed  changes  of  about 
5%  from  Summer  to  Winter  1988-1989.  Fortunately,  these  changes  in 
period  were  slow  enough  that  the  effect  could  be  minimized  by  applying 
a  numerical  correction  to  the  recorded  outputs,  after  computer 
integration. 

(3)  The  output  of  the  torsion  balance  exhibits  a  DC  drift  from 
-10  V  to  +  10V  (a  full  cycle  in  about  two  months).  The  required 
signal,  characterized  by  a  very  small  amplitude  +  (10  microVolt-to-100 
microvolt),'  is  on  top  of  this  slowly  varying  drift.  By  using  a  16  bit 
resolution  in  the  digitization  process  (instead  of  the  12  bit 
resolution  of  Prof.  Weber's  earlier  tests),  we  were  able  to  detect  the 
signal  most  of  the  time. 
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Once  the  relevance  of  these  torsion  balance  difficulties  was  fully 
appreciated,  Raytheon  (with  the  concurrence  of  Prof.  Joe  Weber) 
submitted  a  recommendation  to  DARPA  to  extend  the  present  phase  of  the 
project  for  three-to-six  months.  This  would  at  least,  provide  a  new 
damper  (characterized  by  critical,  or  near-critical,  damping)  and  a 
thermal  control  subsystem  for  the  instrumentation.  Then  the  data 
collection  could  be  resumed  with  26  source-replicas  mounted  on  the  top 
of  the  1  RPM  table.  LANL  is  already  proceeding,  under  direct  DARPA 
funding,  to  fabricate  these  26  replicas  and  to  mount  them  around  the 
rim  of  the  1  RPM  rotating  table,  in  order  to  maximize  the  ratio 
"repulsion  force/attractive  force".  The  repulsion  force  represents 
the  desired  signal,  and  the  attractive  force,  due  to  gravity-gradient 
background,  the  undesired  one.  Computer  simulations,  carried  out  at 
SAO  under  a  Raytheon  subcontract,  have  shown  that,  with  this 
arrangement,  we  could  expect  an  improvement  in  the  ratio  "repulsive 
force/attractive  force"  of  500-1000. 

Thus  far,  data  have  been  collected  by  Raytheon  without  using  the  26 
replicas.  With  all  applicable  caveats,  we  could  draw  the  following 
tentative  conclusions  (we  cannot  stress  strongly  enough  how  tentative 
our  current  results  are): 

(1)  There  is  a  difference  between  the  recordings  (obtained  with 
168  hours  integration  time)  performed  with  tritium  and  with  deuterium. 
This  difference  would  be  consistent  with  a  repulsion  force  which  is 
present  when  experimenting  with  the  tritium.  The  intensity  of  this 
repulsion  force  appears  to  be  of  a  few  microdynes.  This  tracks 
closely  Prof.  Weber's  observations  at  U.  of  Maryland  in  1986; 

(2)  If  a  1/4"  lead  shield  is  wrapped  around  the  8"  diameter 
cylinder  that  houses  the  torsion  balance,  this  repulsion  force  seems 
to  disappear; 
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( 3 )  Observations  ( 1 )  and  ( 2 )  above  have  been  performed  only 
once,  and  the  SNR  was  rather  limited.  That  is  why  we  say  that  the 
experimental  evidence  is,  thus  far,  most  tentative.  Should  these 
results  be  confirmed  by  the  forthcoming  tests  during  the  contract 
extension,  the  question  will  arise  as  to  the  origin  of  a  radiation 
pressure  that  produces  a  torsion  balance  force  that  disappears  when  a 
lead  shield  is  interposed  between  source  and  sensor.  We  do  not  expect 
that  this  question  will  be  easy  to  answer.  A  first  candidate 
radiation  that  comes  to  mind  is  X-ray  radiation.  However,  where  do 
these  X-rays  come  from?  They  cannot  come  from  the  inside  of  the 
tritium  container:  the  thickness  of  its  wall  is  too  large  to  allow 
even  the  smallest  leak.  Another  possible  mechanism  has  been  proposed 
that  would  generate  14.1  MeV  neutrons,  that,  in  turn,  could  produce 
secondary  radiation  affecting  the  torsion  balance.  Clearly,  we  need 
to  add  new  sensors  to  our  instrumentation  at  LANL  such  as  a 
proportional  counter,  a  X-ray  detector,  a  gamma-ray  detector,  a 
neutron  spectrometer,  etc.  Dependent  upon  the  results  of  the 
measurements  above,  we  could  interpose  shields  between  source  and 
torsion  balance,  such  as:  (a)  IR  shield  made  of  several  layers  of 
gold-coated  mylar,  each  with  95%  reflection  coefficient  in  the  IR 
band,  or  an  equivalent  shield;  (b)  Magnetic  shield  (several  layers  of 
flexible  Permag  metglass);  (c)  1/4"  lead  shield,  against  X-rays;  (d) 
Cadmium/boron/ lead  shield  against  14  MeV  neutrons. 


In  addition  to  the  experimental  activity  that  Raytheon  has  carried  out 
at  LANL  with  the  torsion  balance,  we  have  constructed  and  tested  in 
Portsmouth,  RI  (without  using  tritium  sources  of  neutrinos,  thus  far) 
a  cryogenic  force  sensor  that  has  exhibited  a  sensitivity  of  about 
10- *  dyne  in  104  seconds  integration  time.  Testing  this  cryogenic 
instrument  (that  operates  at  4°K,  is  tuned  at  100  Hz,  and  is  suspended 
from  a  tripod  equipped  with  a  2 -stage  vibration  isolator)  at 
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LANL  requires  a  2000  RPM  rotating  table,  with  three  tritium  sources 
mounted  on  the  table  top,  120°  apart,  to  modulate  the  neutrino  flux 
amplitude  at  100  Hz.  This  table  has  not  yet  been  built,  although  an 
analysis  of  its  feasibility  and  a  preliminary  design  have  been 
completed  by  Raytheon  and  SAO. 

Should  the  tests  with  the  torsion  balance,  to  be  completed  by  December 
1989,  provide  support  to  Prof.  Weber's  expectation  of  an  abnormally 
large  coherent  scattering  cross-section  using  tritium- radiated 
neutrinos,  we  would  then  have  a  simpler  and  cheaper  alternative  for 
the  100  Hz  modulation "  of  neutrino  flux.  This  would  consist  of  a 
modified  Weber  chopper,  rotating  at  100/n  RPS,  where  n  is  the  number 
of  sapphire  crystals  mounted  at  the  periphery  of  the  rotating  disk. 
The  chopper  could  be  used  at  LANL  with  a  tritium  source  kept  at  a 
fixed  station,  on  a  stand.  By  replacing  the  tritium  source  with  a 
deuterium- filled  container,  we  could  verify  whether  a  tritium- induced 
signal  would  continue  to  be  present,  or  would  disappear,  when  using 
the  deuterium. 

The  cryogenic  force  sensor  that  Raytheon  has  constructed  provides  a 
sensitivity  improvement  of  more  than  4  orders  of  magnitude  with 
respect  to  the  room- temperature  torsion  balance,  thus  reaching  a 
threshold  of  less  than  10* 1 0  dynes,  with  a  168-hour  integration  time. 
The  instrument  has  the  potential  of  achieving  a  sensitivity  of  the 
order  of  a  picodyne  (10- 12  dynes),  by  adding  cryostat  dilution 
refrigeration  provisions,  which  would  bring  the  operating  temperature 
down  to  4  milli°K. 

*  *  * 

Finally,  concerning  the  laboratory  tests  on  the  feasibility  of  the 
magnetic  interaction  sensor,  Raytheon  has  constructed  a  SQUID 
magnetometer  under  the  1988  Independent  Development 
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Program  (IDP),  and  has  performed  the  following  measurements ,  as 
contractually  required: 

(1)  Determination  of  the  SQUID  system  noise  (this  is  called 
intrinsic  noise,  when  the  SQUID  is  unloaded):  for  frequencies  larger 
than  10  Hz,  our  SQUID  has  an  intrinsic  noise  of  about  10* 6  <£>  Hz- 1 ' 2 
(whereto  =  2.07  10- 1  s  Wb).  This  is  a  factor  of  10  better  than 
commercially  available  SQUIDs.  However,  we  need  a  noise  level  less 
than,  or  equal  to,  10* 1  Hz-1/2  ,  when  the  SQUID  is  loaded  by  the 
target.  Therefore,  we  have  still  a  long  way  to  go.  JASON  estimated 
that  a  loaded  SQUID  would  have,  at  best,  a  system  noise  of  10*  *Jk 
Hz- 1 ' 2  .  - 

(2)  Determination  of  the  maximum  achievable  relative 
permeability  for  the  interaction  target.  We  have  measured  a  jjx  =  106  , 
when  operating  on  the  steep  side  of  the  hysteresis  cycle.  We  require 
fix  =  10*  .  We  have  still  a  factor  of  100  to  gain.  JASON  had  estimated 
that  at  best  we  could  achieve  a  jix  =  103  . 

(3)  Determination  of  the  maximum  achievable  collecting  area.  We 
need  10*  cm2  .  Our  SQUID  stopped  operating  with  a  15  cm2  area.  JASON 
had  estimated  that  at  best  we  could  achieve  1  cm2  .  Here  also,  we  have 
a  factor  of  100  to  gain.  The  use  of  a  superconducting  transformer 
offers  great  promise. 

Our  estimate  of  the  effort  that  is  required  to  complete  the 
preliminary  laboratory  tests  on  the  three  basic  JASON  questions 
concerning  magnetic  sensor  feasibility  is  as  follows:  (a)  Laboratory 
work  on  SQUID  improvement  and  on  target  testing,  still  using  a  1  Kg 
mass  for  the  high-permeability  interaction  target,  will  require  about 
3  man-years  effort  over  a  1  year  period.  (This  work  does  not  require 
tritium  sources  and  can  be  carried  out  at  Raytheon,  Portsmouth,  RI); 
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(b)  Completion  of  the  SQUID  improvement  program,  and  laboratory  work 
with  a  25  Kg  interaction  target,  will  require  about  4  man-years  effort 
over  a  1  year  period  (work  still  to  be  done  in  Portsmouth,  RI). 

Should  the  answers  to  the  three  JASON  questions  provide  support  for 
sensor  feasibility,  we  could  then  move  toward  a  program  phase  that 
includes  construction  of  the  actual  sensor.  'I’his  would  be  used  at 
LANL/TA-33  with  tritium  sources  of  neutrinos  (this  instrument  would 
use  a  250  Kg  interaction  target,  housed  in  a  custom-made  4°K  cryostat) 
for  the  conduct  of  fully-probative  detection  experiments. 

Concerning  an  assessment  of  the  likelihood  of  success,  we  must 
distinguish  between  two  goals .  Achieving  the  three  goals  in  the 
instrumentation  subsystem  performance  identified  by  JASON,  could  be 
characterized,  in  Raytheon's  opinion,  as  better  than  50%,  while  the 
probability  of  detecting  neutrinos  with  the  magnetic  sensor  cannot  yet 
be  assessed  at  this  time.  Notwithstanding  this  uncertainty,  our 
strong  recommendation  is  that  support  be  provided  by  DARPA  to  build 
up,  item  by  item,  the  technology  basis  that  is  a  pre-requisite  to 
fully  develop  this  novel  detector  of  low-energy  neutrinos. 
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EXECUTIVE  SUMMARY 


In  tne  t i rst  two  months  of  program  performance,  Raytheon's  effort  concentrated  in 
planning  ♦he  added-scope  activity,  which  will  consist  of  constructing  instruments 
and  performing  detection  tests  using  a  tritium  source.  The  tritium  is  expected  to 
De  provided  as  a  government  furnished  unit  by  Lawerence  Livermore  National 
Laooratory  (LLNL).  An  initial  draft  of  the  program  plan  has  been  reviewed  by  DARPA 
ana  the  JASON  Study  Group.  Raytheon  expects  to  receive  instructions  from  DARPA  as 
to  redirections  of  activity  priorities,  as  formulated  in  the  above  reviews.  At 
that  point,  the  final  version  of  the  plan  will  be  issued,  and  the  project  activity 
will  begin  in  earnest. 

APPENDIX  A,  attached  to  this  report,  reproduces  the  plan  as  submitted  to  DARPA 
and  JASON.  We  expect  that  any  changes  and  redirection  will  affect  only  a  minor 
portion  of  the  plan,  and  will  be  limited  to  the  calibration  of  the  sensors  used 
in  the  forthcoming  experiments,  and  the  modalities  of  the  data  collection.  To 
compensate  for  the  increase  in  scope  of  some  of  the  tasks,  there  might  be  a 
decrease  in  the  amount  of  work  done  for  the  magnetic  interaction  sensor.  This 
would  also  be  in  deference  to  the  concerns  expressed  in  JASON  Report  JSR-84-105, 
entitlea  "Neutrino  Detection  Primer",  dated  June  1987,  that  has  been  critical 
of  Raytheon's  magnetic  approach  to  neutrino  detection. 

At  this  point,  our  expectations  about  the  content  of  the  final  program  plan  are 
that  most  of  the  effort  will  be  spent  in  constructing  a  high-sensitivity  cryogenic 
force  sensor  (essentially  a  gravity  gradiometer  with  one  of  the  proof  masses 
replacec  by  sapphire  crystals,  operated  at  4°K),  and  a  replica  of  Prof.  Weber's 
torsion  balance  and  tuning-fork  instruments  (both  working  at  room  temperature ) ; 
in  collecting  data  with  the  LLNL  tritium  source  of  neutrinos  (a  3.5"  sphere  with 
a  100  Kiiocurie  source  intensity);  and  in  reducing,  processing  and  analyzing  the 
data.  The  source  w i I  I  be  rotated  on  a  revolving  round-table,  capable  of  variable 


angular  velocities  as  high  as  6000  RPM .  We  also  expect  that  the  final  program 
plan  will  include  the  construction  of  a  scaled-down  magnetic  interaction  sensor, 
with  a  target  of  only  25  Kg,  instead  of  the  f  u  I  I  -sea  I e  target  of  250  Kg.  It  will 
use  a  standard-size  4°K  cryostat,  instead  of  the  custom-made  large-size  cryostat 
tnat  was  designed  as  a  part  of  the  previous  contract.  The  scaled-down  magnetic 
sensor  would  be  used  in  laboratory  tests  of  such  fundamental  performance  parameters 
as  the  maximum  achievable  relative  permeability  of  the  high-mu  target  material, 
and  the  minimum  achievable  system  noise,  when  the  squid  pick-up  loop  is  loaded 
oy  the  target  itself.  These  laboratory  tests  would  directly  address  the  concerns 
expressed  by  JASON  of  the  magnetic  sensor's  ability  to  detect  neutrinos,  and 
would  resolve  the  disagreement  between  JASON/Raytheon  on  the  topic  of  the 
required  integration  time:  JASON  states  that  10^4  seconds  are  necessary,  Raytheon 
that  104  seconds  are  enough.  A  ratio  10^0  between  the  two  estimates  translates 
in  a  ratio  of  1010  between  the  estimates  in  Signa I -to-Noi se  ratio  expected  to 
oe  achieved  in  the  proposed  experiment.  Specific  items  of  disagreement  are: 

-  relative  permeability  of  magnetic  material 
JASON  103 

Raytheon  10®  (by  the  Bozorth  treatment  consisting  of  applying 
thermal  processes  and  mechanical  stresses) 

-  SQUID  pick-up  coil  area 
JASON  1  cm2 

Raytheon  10^  cm2  ( y  the  use  of  a  superconducting  transformer) 

-  SQUID  noise  density 


i  v 


JASON  10"4  fjfT 
Raytheon  10”®  f.'F 


It  can  be  seen  from  the  above  that  these  three  areas,  by  themselves,  account 
for  the  10'0  discrepancy  (10^  x  10^  x  10^  =  TO'®).  The  laboratory  tests  with 
the  scaled-down  instrument  would  provide  values  for  each  of  the  above  parameters 
and  would  serve  as  a  realistic  estimate  for  the  performance  of  the  final, 
full-scale  instrument. 
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APPENDIX  A 


THE  PROGRAM  PLAN  FOR  THE  ADDED  SCOPE  EFFORT,  AS  PRESENTED 
TO  JASON  ON  7/10/1987,  AT  LA  JOLLA,  CA. 


This  Appendix  contains  the  copy  of  the  viewgraphs  used 

•• 

by  Raytheon  in  the  presentation  to  JASON. 
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WEBER  TORSION  BALANCE  -  ELECTRONICS 


SotMATIC  DIAGRAM  OF  OVERALL  INSTRUMENT 


VEBEH  TORSION  BALANCE  -  ELECTRONICS  (continued) 


Schematic  diagram  of  phase-sensitive  detector 
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Integration  time  /v] 0?  sec 
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EXECUTIVE  SUMMARY 


In  che  past  quarter,  activity  for  most  of  the  project  proceeded  close  to  the 
pre-established  schedule,  although  we  are  experiencing  delays  in  some  areas. 

Progress  was  substantial  for  the  cryogenic  force  sensor.  The  construction  of 
the  room-temperature  gravity  gradiometer  at  IFSI-CNR  (to  be  used,  after  modifi¬ 
cation,  by  SAO  and  Raytheon)  was  completed  in  early  October  1987.  The  PMO 
witnessed  a  laboratory  demonstration  of  its  functionality  on  October  13,  1987. 
Sapphire  crystals  (being  procured  from  Union  Carbide  by  Raytheon)  will  be  embedded 
in  one  of  the  two  proof -masses  at  IFSI-CNR  prior  to  shipment  of  the  modilied 
instrument  to  SAO.  This  shipment  is  scheduled  to  take  place  on  or  before  A/30/88. 
In  the  meantime,  SAO  is  scheduled  to  procure  the  4eK  cryostat,  the  tripod  for  the 
support  of  the  instrument,  the  related  suspension,  etc.,  so  that  a  fully  inte¬ 
grated  measurement  system  will  be  made  available  to  Raytheon  by  5/31/88.  This 
will  enable  initiation  of  the  feasibility  tests  with  a  neutrino  source  of  100 
Kilocurie,  in  an  area  suitably  licensed  by  the  Nuclear  Regulatory  Commission  (NRC). 
These  tests  are  expected  to  last  four  months,  until  9/30/88. 

There  was  also  progress  in  the  SAO  design  of  the  rotating  source  assembly,  although 
this  design  could  not  be  completed,  as  planned,  in  this  last  quarter.  LLNL  has 
not  yet  provided  Raytheon  with  drawings  of  the  3.5"  tritium  source  sphere,  in¬ 
cluding  information  on  its  mass.  The  rotating  table,  with  vertical  axis  of 

rotation,  will  turn  at  3000  RPM  and  will  be  characterized  by  s  rotational  stability 
A 

of  1  part  in  10  .  At  3000  RPM,  two  sources,  mounted  at  diametrically  opposite 
points  on  the  rim  of  the  round  table,  will  generate  a  signal  at  100  Hz,  which  is 
the  resonance  frequency  of  the  cryogenic  gravity  gradiometer.  If  three  sources 
were  available,  the  rotational  speed  could  be  2000  RPM,  still  providing  a  signal 
at  100  Hz.  The  reason  for  multiple  sources  is  to  obtain  a  signal  frequency 
different  than  the  frequency  of  mechanical  rotation  of  the  table  (50  Hz  for  3000 
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RPM),  the  latter  being  the  frequency  at  which  the  mechanical  vibration  noise 
produced  by  the  table  is  expected  to  be  at  its  maximum.  An  adaptive  servomechanism 
speed  control  will  make  it  possible  to  choose  other  values  of  rotational  speed, 
such  as  30  RPM  ( 1  Hz  signal  frequency  for  2  sources)  for  experimenting  with  the 
magnetic  interaction  sensor,  or  0.05  RPM  to  0.5  RPM  (1/600  and  1/60  Hz  signal 
frequency  respectively,  again  with  two  sources),  for  experimenting  with  Prof. 
Weber's  torsion  balance.  The  rotating  source  assembly  is  expected  to  be  ready 
for  feasibility  tests  at  Raytheon  on  or  before  29  February  1988,  if  the  missing 
information  from  LLNL  is  available  no  later  than  30. November  1988.  On 
15  March  1988,  we  expect  to  start  experimenting  with  Prof.  Weber's  torsion 
balance  (about  three  months  prior  to  the  availability  of  Raytheon's  cryogenic 
force  sensor). 

In  this  past  quarter,  the  University  of  Maryland  received  from  Raytheon  the 
subcontract  to  construct  a  replica  of  Weber’s  torsion  balance-expected  delivery 
to  Raytheon  on  15  March  1988. 

Concerning  the  tuning  fork  and  the  rotating  shutter,  work  proceeds  on  schedule 
at  University  of  Maryland,  under  a  direct  contract  from  DARPA/AFOSR.  We  also 
understand  that  the  feasibility  tests  with  this  sensor  will  be  carried  out 
directly  by  University  of  Maryland,  without  Raytheon's  involvement. 

Initiation  of  work  for  the  reduced-scope  magnetic  interaction  sensor,  to  be  used 
in  the  laboratory  measurements  of  fundamental  instrument  parameters,  has  ex¬ 
perienced  a  delay  of  about  one  month,  with  respect  to  the  schedule  worked  out  at 
the  planning  meeting  at  IESS-CNR  on  12  October  1987.  The  completion  of  the 
design  of  the  25  Kg  high-permeability  interaction  target  is  now  expected 
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31  December  1987.  The  impact  of  this  delay  on  the  overall  project  schedule  is 
minor,  because  the  laboratory  tests  on  the  magnetic  sensor  represent  an  activity 
that  is  outside  the  main  focus  of  our  present  work.,  which  is  the  verification  of 
Prof.  Weber's  findings. 

At  the  time  of  this  writing,  two  issues  are  of  some  concern  to  the  project: 

(a)  the  lack  of  information,  thus  far,  on  the  LLNL  tritium  source  (this  prevents 
the  completion  of  SAO's  design  work  for  the  rotaclng  source  assembly); 

(b)  the  unresolved  issue  of  the  location  of  feasibility  tests.  Although,  at  this 
point,  site  selection  does  not  hold  back  progress  in  other  areas  of  project 
activity,  it  is  a  necessary  input  for  the  detailed  planning  of  data  collection. 
Since  the  data  collection  is  expected  to  begin  in  March  1988,  the  effort  with  the 
NRC  to  upgrade  a  lab  facility  presently  cleared  to  handle  radioactive  sources 
will  have  to  begin  3  months  earlier  -  no  later  than  15  December  1987.  The  site 
selection  is  a  Raytheon  responsibility,  and  Raytheon  will  concentrate  on  it  in 
the  weeks  ahead.  At  present,  potentially  available  sites  are  Raytheon  Sudbury, 
University  of  Rhode  Island,  and  LLNL  (the  latter  is  a  solution  of  last  resort, 
considering  the  distances  involved). 
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1.  INTRODUCTION 


The  project  has  entered  a  new  and  more  pressing  phase  of  activity:  the 
construction  of  the  instrumentation  Tor  conducting  the  feasibility  tests.  In 
performing  the  design  of  this  instrumentation,  Raytheon  is  taking  into  account, 
under  DARPA  direction,  the  recommendations  formulated  by  JASON  at  the  project 
review  meeting  of  July  10,  1987.  The  consequence  is  that  the  original  design 
must  be  modified  in  several  areas,  with  a  substantial  upgrading  in  the  quality 
of  the  hardware  and  the  rigor  of  the  experimental  methodology  to  be  adopted. 

The  most  significant  changes  are  introduced  in  the  design  of  the  rotating 
source  assembly  (also  called  the  rotating  table,  or  the  rotating  wheel): 

(a)  housing  the  rotating  table  in  a  vacuum  container,  to  prevent  airborne 
mechanical  vibrations  from  reaching  the  sensors; 

(b)  choosing  number  of  sources  N  (to  be  mounted  on  the  rim  of  the  rotating  table) 

such  that  the  frequency  (Hz)  of  the  signal  (  RFM  x  N)  falls  on  a  local  minimum  of 

60 

the  spectrum  of  the  mechanical  vibration  noise  generated  at  the  sensors'  input  by 
the  rotating  table  (this  noise  is  expected  to  have  a  maximum  at  the  fundamental 
rrequency  (RPM/60); 

(c)  reducing  to  the  lowest  possible  value  the  time-variable  gravity  gradient 
forces  produced  by  the  rotation  of  the  table,  and  determine  the  residuals  by  re¬ 
placing  the  N  sources  mounted  on  the  rim  of  the  wheel  with  empty  dummies  (utilizing 
3b  dummies,  instead  of  (36-N)  dummies  and  N  sources); 

(d)  replacing  the  N  sources  in  the  rotating  table  with  dummies  containing  heat 
sources  to  verify  the  extent  of  the  sensor's  response  to  thermal  effectB  (the 
cryogenic  force  sensor  will  be  Immune  from  this  effect,  because  of  the  thick  wall 
of  the  cryostat.  The  torsion  balance,  on  the  contrary,  might  be  affected); 
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(e)  as  above,  but  with  the  dummies  containing  a  permanent  magnet,  to  verify 
the  extent  of  the  sensor's  response  to  tlme-variabLe  magnetic  fields  that  are 
at  the  rotational  signal  frequency.  The  changes  in  the  design  of  the  cryogenic 
force  sensor  are  more  limited  and  consist  of  the  following: 

(a)  embedding,  into  the  aluminum  block  of  one  of  the  so-called  " proof -mas ses " , 
sapphire  crystals  (to  prevent  X-ray  induced  electrostatic  charging  these  crystals 
are  totally  encapsulated  in  the  aluminum);  the  crystals  will  be  in  the  form  of 
ten  rods,  with  a  1  cm  diameter  and  5.2  cm  length.  They  will  be  inserted  in  ten 
cylindrical  wells  drilled  in  the  aluminum  block  of  the  proof -mass. 

(b)  modifying  the  feedback  circuit  of  the  null-servo  that  keeps  the  proof-mass 
at  equilibrium,  to  make  it  possible  to  apply  an  external  electrical  signal  for 
sensor  calibration; 

(c)  providing  a  signal  for  fine  tuning  the  resonance  frequency  of  the  cryogenic 
force  sensor.  This  will  be  derived  from  the  measurement  of  the  rotational  speed 
of  the  rotating  table. 

The  changes  in  the  magnetic  interaction  sensor  represent  a  reduction  in  scope  of 
the  work  planned  for  this  instrument.  This  work  now  consists  of  laboratory 
measurements  or  fundamental  sensor  parameters,  exclusive  of  attempts  to  detect 
magnetization  induced  in  a  high-permeability  target  by  low-energy  neutrinos.  We 
have  already  started  the  design  of  a  25  Kg  interaction  target  (instead  of  the 
contemplated  250  Kg  version)  and  we  will  measure  the  following  system  parameters 
with  a  high-sensitivity  SQUID  (provided  to  SAO  by  IESS-CNR,  Roma,  Italy): 

(i)  thermal  noise  of  the  SQUID  when  its  input  is  loaded,  via  a  superconducting 
transformer,  by  the  25  Kg  target; 
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(2)  maximum  achievable  relative  permeability  of  the  target  material,  when 
operating  on  the  steep  ascending  side  of  the  hysteresis  cycle; 

(3)  maximum  achievable  target  cross-sectional  area. 

In  FY  88,  analytical  work  has  been  limited  to  the  computer  simulation  of  the  ex¬ 
pected  performance  of  the  instrumentation  being  designed.  At  the  JASON  review  of 
7/10/87,  Raytheon  illustrated  its  plans  to  analyze  the  expected  response  of  the 
cryogenic  gravity  gradiometer  to  forces  expected  to  be  generated  by  the  rotation 
of  the  table.  This  work  has  started,  and  a  first  account  of  what  we  have  done 
ie  provided  in  Section  2.1  of  this  Quarterly  Report:  the  expected  neutrino-induced 
signal  shape  is  illustrated  for  different  numbers  of  neutrino  sources  (1  through 
7).  Next  step  wil  consist  of  computing  the  forces  produced  by  gravity  gradients 
due  to  realistically-estimated  non-idealities  in  the  rotating  table's  cylindrical 
symmetry.  For  the  magnetic  interaction  sensor,  an  analysis  of  the  response  of 
this  instrument  to  neutrons  was  performed.  Section  2.2  illustrates  our  findings 
on  this  topic.  This  analysis  is  useful  in  understanding  the  response  of  the  sensor 
to  neutrinos  and  in  exploring  the  possibility  of  a  sensitivity  calibration  of  the 
instrument  by  the  interaction  of  particles  with  the  target. 

The  progress  achieved  in  this  past  quarter  in  instrumentation  design  and  construction 
is  illustrated  in  Section  3,  which  documents  what  was  accomplished  in  the  design  of 
the  rotating  source  assembly,  and  provides  the  status  on  the  development  of  the 
room-temperature  gravity  gradiometer  at  IFSI-CNR,  Frascati,  Italy  (instrument 
functionality  was  demonstrated  to  Raytheon  on  13  October  1987). 
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2.0  ANALYTICAL  DEVELOPMENTS 


2.1  Signal  Shape  For  Different  Numbers  Of  Neutrino  Sources  On  A  Rotating 
Wheel(*) 


Figure  1.  Neutrino  Sources  on  a  Rotating  Wheel 


Figure  1  shows  a  set  of  neutrino  sources  equally  spaced  around  a  circle 
of  radius  r.  There  is  a  detector  located  at  a  distance  l  from  the  center  of  the 
wheel  on  the  r  -  axis  which  measures  the  i  —  component  of  the  neutrino  flux  from 
the  n  sources  on  the  wheel.  We  wish  to  calculate  the  shape  of  the  signal 

'Contributed  by  David  A.  Arnold,  SAO 


measured  by  the  detector  as  a  function  of  the  rotation  angle  of  the  wheel. 
The  positions  of  the  neutrino  sources  as  a  function  of  time  are  given  by  the 
equations 


where 


and 


Xi  as  r  cos  9% 
W  =  r  sin  6{ 


8%  —  (•'  ~  1)  A8  ■+•  uxt 


A8  —  2r /n 


(1) 

(2) 

(3) 

(4) 


The  source  on  the  z  axis  closest  to  the  detector  is  labelled  t  =  1,  and  the  sources 
are  numbered  counterclockwise.  The  wheel  rotates  with  angular  velocity  u.  The 
distance  of  the  i'**  source  from  the  detector  is 


di  =  \J  Axi  +  |K  (5) 

where 

Axi  *  l-xi  (6) 

The  intensity  of  the  neutrino  flux  at  the  detector  from  the  i>tA  source  is 

/,•  *  C  /£  (7) 


where  C  is  a  constant  and  di  is  the  distance.  The  sum  of  the  z  -  component » 
of  the  fluxes  from  all  the  sources  is 


(8) 


F,  =  V  F, 

« 

A  computer  program  called  FLYWHEEL  has  been  written  to  calculate  the 
value  of  Ft  vs.  the  rotation  angle  of  the  wheel.  A  set  of  runs  has  been  for 
s  =  1  to  7  with  Ft  plotted  at  3  degree  intervals,  t»mg  C  =  100.,  r  =  45  cm  and 
t  *  70  cm.  The  value  of  C  was  chosen  arbitrarily  to  give  numbers  close  to  unity. 
The  flywheel  is  expected  to  be  one  meter  in  diameter  with  spherical  sources  9  cm 
in  diameter.  The  container  for  the  detector  is  expected  to  be  about  14  inches 
(35.56  cm)  in  diameter. 

Figure  2  shows  the  shape  of  the  signal  for  a  single  neutrino  source.  At  the 
closest  approach  the  distance  is  25  cm  and  the  signal  from  equation  (7)  is  .16000  in 
arbitrary  units.  At  the  furthest  point  from  the  source  (urt  =  180  deg)  the  signal  is 
.00756.  The  difference  is  .15244  units.  At  vt  *  90  deg  the  magnitude  of  the 
signal  (with  d  =  83.216  cm)  is  .01444  and  the  x  -  component  is  .01215.  figure 
2a  shows  the  flux  Fa  vs.  rotation  angle  of  the  wheel  at  3  deg  intervals  from  0  to 
360  deg.  Figure  2b  shows  the  signal  from  -180  to  +180  deg.  Figures  3  to  8  show 
the  signal  for  n  *  2  to  7.  The  pulse  width  for  n  =  1  is  about  40  deg.  at  half 
amplitude.  At  this  width,  several  sources  could  be  used  before  there  would  be 
significant  overlap  of  the  signals.  As  more  sources  are  added  the  signal  becomes 
more  sinusoidal. 

Table  1  lists  the  m»rimiim  signal  the  minimum  signal  -  and  the 
difference  AF  between  the  maximum  and  minimum  signals  for  each  of  the  seven 
cases. 


n 

Fmax 

Table  1 

Faun 

AF 

1 

.16000 

.00756 

.15244 

2 

.16756 

.02429 

.14327 

3 

.17829 

.04852 

.12977 

4 

.19185 

.07897 

.11288 

5 

.20808 

.11359 

.09449 

6 

.22681 

.15033 

.07648 

7 

.24780 

.18762 

.06018 

Table  1.  Amplitude  of  the  signal  vs.  number  of  sources  n. 

Table  2  gives  some  of  the  details  of  the  calculations  of  Fa  for  n  s  6  and 
0.  The  geometry  for  this  case  is  given  in  Figure  1. 


Table  2 


i 

6 

X 

V 

d 

F 

F, 

I* 

1 

0 

45 

0.00 

25.00 

.16000 

.16000 

.16000 

2 

60 

22.5 

38.97 

61.44 

.02649 

.02048 

.18048 

3 

120 

-22.5 

38.97 

100.37 

.00993 

.00915 

.18963 

4 

180 

-45. 

0. 

115.00 

.00756 

.00756 

.19719 

5 

240 

•22.5 

•38.97 

100.37 

.00993 

.00915 

.20633 

6 

300 

+22.5 

-38.97 

61.44 

.02649 

.02048 

.22681 

Table  2.  Calculation  of  the  signal  for  n  =  6  and  u it  *  0. 


The  neutrino  sources  to  be  used  In  this  experiment  have  a  strength  of  100  Kilocurie. 

Since  1  Curie  translates  into  3.4  10^  neutrinos/ sec,  we  hav*  an  intensity  of  3.4  10^ 

neutrinos/sec  with  each  source.  The  constant  C  in  equation  (7)  can  be  evaluated  from 

the  equation  ,  .  in15  u 

C  ,  *  2.7  10 

4  y 

In  the  calculations,  a  value  of  100  was  used  to  produce  numbers  on  the  order  of  unity. 

The  flux  values  can  be  converted  to  absolute  units  for  this  source  strength  by  multi- 
12 

plying  by  2.7  10  .  For  example,  with  one  source  at  a  distance  of  25  cm,  the  flux 

11  -1-2 

in  arbitrary  units  was  0.16  and  the  absolute  value  is  4.3  10  neutrinos  sec  cm  . 

The  momentun  per  second  passing  through  the  target  Is  given  by  the  formula 

E, 

j  A 

c 

-1  -2 

where  j  is  the  flux  in  neutrinos  sec  cm  ,  Ey  Is  the  average  energy  per  particle,  c 

is  the  velocity  of  light  and  A  is  the  area  of  the  target.  If  the  average  energy  of  the 

neutrinos  is  10  keV,  the  energy  per  particle  is  1.6  10“8  ergs  (  with  1  eV  ■  1.6  10'12 

2 

ergs).  Using  a  target  area  of  52  cm  gives 

j  EV  A  »  1.2  10”8  dynes, 
c 

This  is  the  force  on  the  target  if  all  the  momentum  of  the  incident  beam  is  transferred 
to  the  target.  In  the  arbitrary  units  used  in  the  signal  strength  calculations,  the 
signal  was  0.16  for  the  case  of  one  source  at  25  cm.  The  arbitrary  units  can  be  convert¬ 
ed  to  dynes  by  multiplying  by 

1.2  10'5  , 

-  ■  7.5  10  . 

0.16 

The  neutrino  flux  is  Incident  upon  a  target  which  Is  a  high-Q  mechanical  oscillator. 

The  behaviour  of  the  detector  can  be  computed  from  the  mass,  restoring  stiffnes,  and  damp¬ 
ing  of  the  target.  The  force  as  a  function  of  time  is  given  by  curves  such  as  Figures 
2  through  8.  A  computer  program  has  been  written  to  calculate  the  response  of  the  detector 
for  various  input  signal  shapes.  Calculations  are  also  being  done  to  evaluate  the  gravity 
force  on  the  detector  from  the  rotating  masses.  The  results  of  these  studies  will  be 
presented  in  the  next  quarterly  report. 
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3.  INSTRUMENTATION  HARDWARE  EFFORTS 


3.1  Conceptual  Design  of  the  Rotating  Source  Assembly  ^ 

The  rotating  source  assembly  Is  a  circular  table,  having  a  diameter  of  about  1  meter, 
that  revolves  around  a  vertical  axis  at  a  speed  variable  from  0.05  RPM  to  3000  RPM.The  speec 
is  adjustable  By  external  control  and  Is  stabilized  within  1  part  in  10*  by  an  active  servo- 
loop  mechanism.  Figure  9  shows  that  the  rim  of  the  table  accomodates  the  neutrino  sour¬ 
ces,  in  number  that  can  be  chosen  from  2  to  4.  The  source  is  a  sphere  with  about  9  cm 

diameter.  If  there  are  two  sources  in  the  rim,  they  are  diametrically  opposite  and  the 

space  in  between  is  filled  with  Identical  spheres  that,  however,  are  dummies  and  do  not 
contain  tritlun  in  their  interior.  We  expect  to  use  a  number  of  dunmies  from  a  minimum  of 
22  to  a  maximum  of  34  (  assuming  to  use  two  sources,  which  makes  for  a  total  number  of 

spheres  in  the  rim  from  24  to  36).  The  objective  that  we  expect  to  achieve  with  the  arrange 

ment  above  is  as  follows  : 

(a)  We  can  have,  as  signal  frequency,  a  value  as  high  as  100  Hz  (  for  instance 

rotating  two  sources  at  3000  RPM,  or  three  sources  at  2000  RPM).  By  doing  so,  we  simplify 

tremendously  the  protection  of  the  sensor  from  vibration  noise,  because  a  rather  simple 

2 

suspension  that  resonate  at  1  Hz,  can  Immediately  provide  an  attenuation  of  0 00/1 1  »or 
80  dB  Ctwo  in  series  will  give  160  dB  attenuation).  We  would  lose  totally  this  attenua¬ 
tion,  if  the  signal  frequency  would  be  1  Hz  instead.  This  is  an  adavantage  that  can  be 
hardly  overemphasized,  as  it  is  well  known  to  every  scientist  who  works  with  force  sensors 
at  frequency  in  the  vicinity  of,  or  below,  1  Hz; 

(b)  By  filling  the  rim  with  dunmies  identical  in  all  respects  to  the  source,  but 
devoid  of  tritium,  we  minimize  the  time-variable  gravity  gradient  forces  produced  by  the 
rotation  of  the  table.  The  residual  time-variable  component  will  have  a  spectrum  centered 
at  a  frequency  — —  x(Nunber  of  dunmies  +  Number  of  sources),  therefore  well  outside 
the  passband  of  the  receiver  (  tuned  at  100  Hz,  and  with  a  bandwidth  of  10~*  Hz). 

(c)  The  rotating  table,  with  the  rim  suitable  for  accomodating  sources  and  dunmies 

can  be  utilized  to  carry  out  measurement  of  the  sensitivity  of  the  instrunentation  under 

test  to  heat  sources  and  to  magnetic  sources. This  can  be  done  by  placing  into  a  dunmy  a 

heat  generator  (  for  instance  a  resistor  fed  by  a  battery)  or  a  permanent  magnet.  By 

rotating  the  table,  we  can  vary  the  frequency  of  these  spurious  signals  and  check  the 
response  to  them  of  each  sensor. 


(*)  Contributed  by  Dr.  M.D.Grossi,  P.I.,  Raytheon  Co. 
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3.2  A  first  example  of  a  possible  mechanization  of  the 
rotating  source  assembly  v 


The  requirements  for  the  rotating  source  assembly  can  be  recapitulated  as 
follows  : 

1)  Speed  must  be  variable  from  an  upper  limit  of  3000  RPM  (  if  at  all  possible. 
6000  RPM)  to  a  lower  limit  that  Is  close  to  zero  (  something  like  0.05  RPM); 

2)  We  need  the  ability  to  set  the  speed  of  the  table  to  match  the  resonant 
frequency,  or  the  most  suitable  frequency,  of  one  or  more  detectors  used  in  the  feasibility 
tests.  The  cryogenic  force  sensor  (gravity  gradiometer)  has  Its  resonance  at  approximately 
100  Hz,  and  this  determines  the  maximum  frequency  of  operation.  The  resonance  is  fairly 
sharp—  a  Q  of  about  10,000—  so  the  turntable  frequency  must  have  a  stability  of  about 
0.001  Hz  p-p  variation  if  it  Is  not  to  jitter  in  and  out  of  resonance; 

3)  At  least  one  of  the  detectors  requires  a  frequency  reference  signal  for  lock¬ 
ing  its  internal  wheel  to  the  source  rotation  frequency; 

A  system  to  meet  all  these  requirements  can  be  shown  in  Figure  10.  In  the  Figure, 
the  turntable  is  phase-locked  to  a  precision  variable  frequency  generator,  which  also  drives 
a  frequency  synthesizer,  which  in  turn  provides  a  reference  signal  for  any  sensor  system 
that  requires  tt. 

fo 

A  literature  search  was  done  among  all  the  product  literature  pertainingAhigh-torque 
variable  speed  motors.  Motors  in  the  6000  RPM  range  seem  to  have  insufficient  torque.  The 
only  off-the-shelf  motor  that  could  be  identified  was  MAVILOR  Model  M010.000,  with  rated 
torque  of  375  lb-in  at  a  speed  of  3000  RPM.  Figures  11a  and  lib  are  catalog  pages  describ¬ 
ing  this  motor.  The  manufacturer  also  provides  a  driver  module  for  the  motor,  capable  of 
supplying  200  V  DC,  at  60  A.  This  driver  has  some  digital  control  capability,  so  it  may  be 
feasible  to  design  In  a  small  single-board  computer  to  handle  all  the  loop  control  and  fre¬ 
quency  compensation. 

When  we  know  more  about  the  nunber  of  sources  that  we  will  utilize  in  the  rim  of  the 
rotating  table,  we  will  be  able  to  formulate  a  more  complete  set  of  design  requirements  and 
do  a  more  thorough  study  of  the  motor  market  before  we  can  make  a  hardware  selection. 


C*1  Contributed  by  Larry  Coyle,  SAO  Central  Engineering. 
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M  AVI  LOR  Motors 

Model  MOI 0,000 


Compact  “flat-style”  construction  permits  Mavilor  Series  10,000  DC  servo¬ 
motor  to  produce  up  to  18  ftp  In  lees  than  In.  of  motor  length.  Rated  for 
continuous  duty  at  3,000  rpm  with  a  maximum  speed  of  3£00  rpm,  they  will 
run  as  slowly  as  0.1  rpm  under  load  without  cogging,  and  with  infinitely 
variable  control  of  all  speeds  in  between. 

The  ironless  dish-shaped  armature  limits  inertia  to  l/10th  that  of  a  con¬ 
ventional  motor  and  holds  inductance  to  less  than  200  /uH.  Low  Inertia  and 
high  instantaneous  torque  make  the  Mavilor  motor  ideal  tor  applications 
that  require  high  acceleration  and  deceleration  and  fast  reversals.  Low  in¬ 
ductance  and  planar  commutation  prevent  arcing,  extending  brush  Ilfs  to 
over  10,000  hr.  Copper  wire  windings  embedded  in  epoxy  form  a  strong, 
rigid  armature  for  long  life  in  demanding  service. 


Specifications 


Model 

MO1O000 

Fan 

1 

8 

Cont.  Torque  Rating 

Tc 

02-In. 

5987 

Ib-in. 

375 

Output  Power' 

Po 

hp 

17.8 

Rated  Voltage  (*5%) 

Vr 

Vdc 

237 

Rated  Current 

Ir 

amps 

81 

Rated  Speed 

Nr 

rpm 

3000 

Max.  Power  Rate 

A 

kW/sec 

2300 

Efficiency 

% 

92 

Pnak  Torque  Rating* 

Tp 

oi-ln. 

48144 

Ib-ln. 

3009 

Maximum  Speed* 

Nm 

rpm 

3500 

Torque  Sensitivity  (*5%) 

Kt 

02  In./ amp 

102 

Volts  Back  EMF  (sS%) 

Ke 

Volts/ krpm 

78 

Friction  Torque 

Tf 

oa-ln. 

42.6 

Viscous  Friction 

Fi 

02-ln./krpm 

56.6 

Terminal  Resistance 

Ra 

ohms 

.14 

Armature  Inductance 

La 

AJM 

196 

Moment  of  Inertia 

Jm 

02-in.-sec* 

7.08 

Mechanical  Time  Constant 

tm 

ms 

20 

Continuous  Service  Characteristics 

Model  MO  lOjOOO 

£>  -«r  ■mjawgv  JT 

xoo  £. 

aooo  aSsigfr  ,Se 


■  sit  j-b 


ToiMIMU  m 


sn  sao 


Notes: 

'  With  a  27.6  x  27.6  x  0.6  inch  metal  heat  sink. 

1  For  maximum  2.0  seconds  at  a  1%  duty  cycle. 

1  Maximum  speed  in  continuous  service. 

At  no  load  can  run  at  4000  rpm  for  5  minutes. 


Figure  11  a 

MAVILOR  Motors,  Model  M010.000 
(Specifications) 


-Z5- 


Dimensions  in  inches  (millimeters) 


Thermal  time  constant,  sec 

Radial  load',  lb 

242 

—armature  to  case 

390 

Axial  load,  lb 

165 

—case  to  ambient 

Thermal  resistance,  *C/Watt 

2250 

Weight,  lb 

187 

—armature  to  case 

.11 

'Radial  load  applied 

-case  to  ambient 

.03 

at  midlength  of  shaft 

achometer  Option 

This  tachometer,  specifically  designed  for  integration 
into  Mavilor  Motors,  has  a  very  low  inertia  and  low  ripple 
assuring  optimum  performance. 


Accessories  and  Options 

To  fulfill  broad  design  requirements,  the  Mavilor  Motors 
can  incorporate  integrally  mounted  and  industrially 
housed  accessories: 

•  Tachometer  •  Encoder 

•  Brake  •  Qearhead 

Other  options  include  dual  shaft,  special  shaft,  hollow 
shaft,  mounted  speed  reducer  and  brush  monitoring. 


Characteristics 

Voltage  constant,  v/krprn  to 

Voltage  ripple,  %  1.5 

Voltage  linearity  at  3600  rpm.  %  0.15 

Bidirectional  tolerance,  %  0.10 

Temperature  coefficient,  V./*C  -0.02 

Inertia,  ox-in.-tec*  3.77x10-* 

Armature  resistance,  ohms  _  98.4 

Inductance,  mH  30 

Current,  mA  5 

Maximum  current,  mA  30 

Maximum  velocity,  rpm  6000 

Armature  weight,  oz  4.8 

Stator  weight,  oz  8.3 


INFRANOR  Drive  Systems 

Infranor  supplies  complete  servo  systems  comprising 
appropriate  Mavilor  motors  with  wideband  PWM  ampli¬ 
fiers  and  other  controls  and  components  Incorporating 
any  available  options  and  accessories  appropnate  to 
specific  applications. 


INFRANORI 

285  Murphy  Road  •  Hartford,  CT  06114 
Telephone:  (203)  525-7743  •  Telex:  96-6469 


Figure  11  b 

MAVILOR  Motors,  Model  MO 10, 000 
(Dimensions  and  options) 

fcvweew  Mint  a 


3.3  Status  Report  on  the  development  of  the  room-temperature 

{*)“  “  . 

gravity  gradiometer  at  IFSI-CNR  '  ' 


On  October  13,  1987,  the  PMO  attended  a  project  review  at  IFSI-CNR,  Frascati, 
Italy  ,  and  witnessed  laboratory  tests  on  the  room- temperature  gravity  gradiometer 
developed  by  Prof.  F.  Fuligni  for  use  at  SAO  and  at  Raytheon  in  the  neutrino  detector 
feasibility  tests.  The  performance  parameters  were  as  follows: 

•  Mass  250  grams  . 

•  Temperature  Room-temperature 

•  Resonance  Freq.  About  22  Hz  without  negative  spring 

"  15  "  with 

•  Q  factor  1.5  10^ 

•  Weakest  7 

detectable  force  4  10”'  dynes  (  without  vibration  abatement^ 
and  with  10“*  conrnon  mode  rejection) 

The  gradiometer  above  is  now  undergoing  modifications  to  become  suitable 

for  use  in  our  program.  We  must  add  sapphire  crystals  to  one  of  the  two  proof -masses 

and  Figure  12  showra  possible  wy  of  doing  this.  The  Figure  illustrates  how  ten  rods 

of  sapphire  could  be  embedded  inside  the  aluminun  block  of  the  proof  mass,  with  no 

parts  of  the  rods  being  exposed  outside  the  aluminun.  This  is  to  prevent  electrostatic 

charging  of  the  rods  due  to  possible  X-ray  Illumination.  Each  rod  is  a  cylinder  with 

5.2  length  and  1  cm  diameter.  The  mounting  of  the  rods  into  the  proof -mass  will  be 

done  in  such  a  way  that  the  Q  of  the  mechanical  oscillator  will  not  be  affected  by 

this  addition.  The  resonance  frequency  will  be  Increased  to  100  Hz. 

At  SAO,  this  gravity  gradiometer  will  be  mounted  inside  a  4°K  cryostat  for  use 

in  the  feasibility  tests.  We  expect  that,  by  cooling  the  instrument  to  liquid  Helium 

4  4 

temperature,  the  Q  will  Improve  and  become  10  ,  so  that,  with  a  10  seconds  Integration 

ft 

time,  we  will  be  able  to  detect  at  threshold  a  force  of  1.3  10  dynes  (  DARPA  require- 

o 

ment  is  to  achieve  the  design  goal  of  10  dynes). 

Figures  13  through  16  show  several  views  of  the  instrument  that  was  tested  at 
IFSI-CNR  on  10/13/87.  Figure  13  depicts  the  gradiometer  baseline,  with  the  two  sensing 
cells  at  Its  ends.  Figure  14  shows  the  electronic  subsystem  with  its  cover  removed;  ft 
contains  the  bridge,  the  feedback  system,  the  preamplifier,  low-pass  filter,  etc.  In 


(*)  Contributed  by  M.O.Grossi,  with  Inputs  from  F.  Fuligni,  TFSI^CNR, 


the  cryogenic  version  of  the  instrument,  to  be  implemented  at  SAC),  the  electronics 
will  be  housed  external  to  the  cryostat.  The  reason  is  that  the  system  noise  is 
determined  only  by  the  temperature  of  the  proof-mass  .that  is  at  4°K  .  In  fact  ,  the 
sensor  is  a  mechanical  oscillator  that  works  at  resonance  :  under  these  conditions, 
the  contribution  of  the  preamplifier  noise  is  negligible.  This  is  also  the  reason  why 
the  utilization  of  quantum  non-demolition  approaches  is  not  required  here  (QND  is 
useful  when  the  preamplifier  noise  predominates  and  we  want  to  avoid  the  back-action 
on  the  oscillating  mass  ).  Figure  15  gives  further  details  of  the  electronics  package, 
while  Figure  16  shows  the  specular  quality  of  the  surface  finishing  by  high-precision 
machine  tools. 


Figure  13  -  The  gradiometer  base 
sensing  celis  at  its 


Figure  14  -  The  electronic  subsystem  with 
Its  cover  removed 


Figure  15  -  A  closer  look  at  the  electro 
nics  package 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  project  has  entered  a  hardware  fabrication  phase,  where  a  significant 
coordination  effort  is  necessary  to  integrate  hardware  from  institutions  as 
IFSI-CNR,  IESS-CNR  (both  in  Italy),  University  of  Maryland,  SAO,  LLNL,  Meyer 
Tool  Mfg.  Co.,  Union  Carbide,  Arnold  Engineering,  etc.  We  feel  that  we  are 
closely  tracking  the  critical  paths,  with  the  expectation  of  soon  resolving  two 
present  causes  of  concern;  the  availability  of  the  tritium  source  drawings  (and 
the  source  itself),  and  selection  of  a  NRC-licensed  test  site.  Of  great 
encouragement  was  witnessing  the  demonstration  of  the  room- temperature  gravity 
gradiometer  in  October,  at  IFSI-CNR,  by  Prof.  Fulignl  and  Dr.  Bordonl  and  to -aee 
the  enthusiasm  and  the  competence  by  Prof.  Carelli  and  his  assistants  at  IESS-CNR. 
They  are  preparing  themselves  to  arrive  at  SAO  with  a  SQUID  that  is  at  least  two 
orders  of  magnitude  better  than  any  commercially  available  SQUID. 

With  some  help  from  DAKPA  to  break  the  administrative  logjam  to  obtain  the  3.5” 
tritium  source,  we  are  confident  chat  we  will  recover  the  current  one-month 
delay  in  this  task. 
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EXECUTIVE  SUMMARY 


The  program  activity  in  the  quarter  reported  herewith  has  been  truly  substantial 
and  the  definition  of  the  feasibility  experiment  has  progressed  significantly. 

Also  visible' progress  has  been  made  in  the  construction  of  the  sensors  to  be  used 
in  the  experiment. 

The  torsion  balance  has  been  constructed  at  the  University  of  Maryland  and  is 
now  undergoing  tests  there,  under  the  direction  of  Prof.  Joe  Weber,  by  him  and 
his  associates.  This  sensor  is  expected  to  be  delivered  to  Raytheon,  Portsmouth, 
Rhode  Island,  sometime  during  April  1988. 

The  cryogenic  force  sensor  is  also  under  construction  at  IFSI-CNR,  Italy  and  the 
sapphire  rods  have  been  already  mounted  in  one  of  the  proof  masses.  This  sensor 
is  expected  to  arrive  in  Portsmouth,  R.I.  in  mid  June  1988,  after  system  tests  at 
I FS I -CNR  inside  a  4°K  cryostat,  to  be  procured  by  IFSI-CNR  in  Europe,  and  to  be 
then  shipped  to  Portsmouth,  R.I. ,  together  with  the  2-cell  force  sensor. 

Concerning  the  tuning  fork,  that  Prof.  Weber  is  constructing  under  a  direct  contract 
from  DARPA,  we  understand  that  it  is  in  advanced  stage  of  construction.  Raytheon  is 
preparing  the  test  area,  ready  to  accommodate  also  the  tuning  fork,  disregard  of  the 
final  disposition  by  DARPA  of  this  sensor. 

Raytheon  is  also  setting-up,  with  in-house  funds,  a  SQUID  magnetometer  to  be  added 
to  the  array  of  sensor  in  the  test  area,  and  to  be  used  to  fulfill  the  assignment 
consisting  of  measurement  three  basic  parameters  of  the  magnetic  interaction  sensor 
(highest  achievable  relative  permeability  in  the  interaction  target,  largest 
achievable  pick-up  area,  and  lowest  achievable  intrinsic  noise  density  when  the 
SQUID  is  loaded  by  the  large-mass,  high-permeability  target). 

Should  DARPA,  in  early  May,  provide  Raytheon  with  a  suitable  tritium  source  (now 
expected  to  be  of  cylindrical  shape),  Raytheon  would  be  able  to  start  experimenting 
in  mid  May  1988,  with  the  following  set-up: 

(a)  Torsion  Balance,  with  the  cylindrical  source  mounted  on  a  very  simple, 

1  RPM  rotating  table,  equipped  with  an  optical  sensor  of  its  angular 
position  and  rotational  speed.  For  this  initial  test,  Raytheon  needs  a 
tritium-filled  cylinder  and  2  dummy  cylinders.  We  understand  that  these 
three  items  would  be  available. 

(b)  Tuning  Fork,  with  the  cylindrical  source  kept  in  a  fixed  position, 
mounted  on  a  pedestal  (replaceable  with  a  dummy,  in  verification  tests), 
to  be  used  in  association  with  a  sapphire-crystal-equipped,  high-speed 
rotating  shutter  (to  be  built  by  Prof.  Weber  under  direct  contract  from 
DARPA  )  .interposed  between  the  tritium  source  and  the  tuning  fork,  in 
order  to  function  as  the  alleged  modulator  of  the  neutrino  beam. 


In  mid  June  1988,  Raytheon  would  add  to  the  two  sensors  in  items  (a)  and  (b)  above, 
the  cryogenic  force  sensor.  This  could  be  tested  just  by  replacing  the  tuning  fork 
with  it,  leaving  everything  else  unchanged. 

In  mid  July,  then,  Raytheon  is  expected  to  receive  from  SAO  the  high-speed  (2000 
RPM)  rotating  table,  with  air  bearing  and  vacuum  enclosure  (subject  to  funding 
authorization  by  DARPA  for  the  added  scope),  to  be  equiped  with  three  spherical 
sources  (  3.5"  spheres,  or  equivalent),  and  with  30  dummies  for  gravity  gradient 
signal  minimization.  The  rotating  table  will  make  it  possible  to  test  the  cryogenic 
force  sensor,  as  well  as  the  tuning  fork,  by  the  variable  neutrino  flux,  due  to 
distance  modulation.  Initial  funding  for  the  construction  of  this  high-speed  rotating 
table  is  about  to  be  awarded  to  SAO  by  Raytheon,  while  the  full  added-scope  contract 
is  still  under  negotiation,  at  this  time. 

Raytheon  plans  to  carry  out  the  feasibility  experiments,  inclusive  of  Weber  verifi¬ 
cation  tests,  in  its  Laboratories  in  Portsmouth,  R.I.  The  specific  laboratory  area  is 
about  to  be  selected  among  three  possibilities.  We  are  carrying  out  measurements,  by 
use  of  accelerometers,  of  the  intensity  of  mechanical  vibrations  on  the  floor  of  each 
candidate  area  and  we  will  select  the  quietest  location.  We  have  also  initiated  dis¬ 
cussions  with  the  Rhode  Island  Nuclear  Regulatory  Agency,  with  the  objective  of  obtain¬ 
ing  from  themthe  license  to  handle  radioactive  materials.  Raytheon  fully  expects  that 
this  license  will  be  obtained  in  time  for  the  mid-May  tests. 

Still  of  some  concern  is  the  fact  that  Raytheon  has  not  yet  received  from  LLNL  the 
needed  information  on  the  size,  mass  and  precise  characteristics  of  the  tritium  sources. 
This  lack  of  information  is  preventing  SAO  from  completing  the  design  of  the  rotating 
table.  We  are  proceeding  here  at  a  risk,  with  making  assumptions  that  might  not  be 
confirmed,  once  that  the  sources  arrive,  and  that  might  require  extra  work. 

On  the  analytical  side,  we  have  completed,  during  the  report  period,  the  analysis  of 
of  tne  resDonse  of  the  cryogenic  force  sensor  to  various  solicitations  produced  by 
the  rotating  table  and  its  distribution  of  sources  and  dummies.  This  quarterly  report 
is  totally  devoted  to  the  illustration  of  the  results  of  this  analysis. 
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1.0  INTRODUCTION 

We  devote  this  quarterly  report  to  the  study  performed  at  SAO  by  David 
A.  Arnold  on  the  subject  of  the  performance  of  the  high-speed  rotating  table  (also 
called  the  rotating  wheel),  and  concerning  the  design  of  the  experiment  to  be 
carried  out  with  the  cryogenic  force  sensor  installed  in  the  vicinity  of  the  wheel, 
near  its  edge. 

In  the  last  quarterly  report,  analyses  were  presented  that  gave  the  shape 
of  the  signal  to  be  expected  from  a  set  of  £  sources  symmetrically  placed  on  the 
rotating  wheel.  The  present  report  presents  the  results  of  analyses  of  the  detector 
to  various  signal  inputs. 

The  detector  has  been  modelled  both  as  a  single  oscillator  and  as  a  compound 
system,  consisting  of  both  a  detector  and  a  case  which  is  free  to  move.  The  motion  of 
the  case  is  particularly  significant  for  gravitational  forces.  Both  analytic  and 
numerical  solutions  have  been  obtained  for  the  one  and  two  mass  models.  Studies  have 
also  been  done  of  the  influence  of  the  wheel  diameter  on  both  neutrino  and  gravity 
signals. 

Several  software  codes  have  been  worked  out  and  have  been  used  in  connection 
with  plotting  subroutines.  The  results  are  displayed  in  Figures  1  through  20.  The 
captions  of  these  figures  are  as  follows: 

Figure  1.  Resonant  response  with  1  percent  of  critical  damping. 

Figure  2.  Resonant  response  with  critical  damping. 

Figure  3.  Resonant  response  with  no  damping. 

Figure  4.  Response  with  the  driving  frequency  1.5  times  the  resonant  fre¬ 

quency. 


Figure  5. 
Figure  6. 
Figure  7. 
Figure  8. 
Figure  9. 
Figure  10. 
Figure  11. 
Figure  12. 
Figure  13. 

Figure  14. 
Figure  15. 
Figure  16. 
Figure  17. 

Figure  18. 


Response  to  a  single  rotating  source  with  a  fast  risetime  detector. 
Resonant  response  to  a  single  rotating  source  with  no  damping. 
Resonant  response  to  a  sinusoidal  input  with  amplitude  .5. 
Resonant  response  with  2  sources. 

Resonant  response  with  3  sources. 

Resonant  response  with  4  sources. 

Resonant  response  with  5  sources. 

Resonant  response  to  a  cosine  driving  force  of  amplitude  .5. 

Response  with  a  driving  frequency  1.5  times  the  resonant  frequency 
and  Zo  —  -  .8. 

Resonant  response  with  10  percent  of  critical  damping. 

Response  with  u  =  .5 14,  and  Zo  =  4/3. 

In  text. 

Resonance  curve  for  mass  2  with  Fi  =  0,  and  F2  =  mj. 

a.  Real  part  of  At  vs.  frequency. 

b.  Imaginary  part  of  At- 

c.  Magnitude  of  At 

d.  Phase  of  At 

Resonance  curve  for  mass  2  with  Fi  =  mi,  and  Ft  =  .86  m2. 

a.  Real  part  of  At  vs.  frequency. 

b.  Imaginary  part  of  At- 

c.  Magnitude  of  At- 

d.  Phase  of  At- 
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Figure  19.  Resonance  curve  for  mass  2  with  F\  -  mi,  and  Ft  =  mt. 

a.  Real  part  of  At  vs.  frequency. 

b.  Imaginary  part  of  At. 

c.  Magnitude  of  At  (638.5  to  639.0  hertz). 

d.  Magnitude  of  At  (637.8  to  639.8  hertz). 

e.  Phase  of  At. 

Figure  20.  Response  of  the  system  with  equal  masses,  equal  driving  forces,  and 
about  10%  of  critical  damping. 

a.  Position  of  mass  1  vs.  time. 

b.  Position  of  mass  2  vs.  time. 
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2.0  RESPONSE  OF  THE  DETECTOR  FOR  VARIOUS  SIGNAL  INPUTS* 

The  neutrino  flux  from  the  rotating  tritium  sources  is  incident  on  a  detector 
which  is  a  high  Q  mechanical  oscillator.  The  frequency  of  the  oscillator  is  tuned 
to  that  of  the  input  signal.  The  motion  of  the  detector  is  described  by  the 
equation 

mz+  bi  +  kz  =  f(t)  (1) 

A  program  called  DETECT  has  been  written  to  solve  equation  (1)  numerically 
using  the  Gear  integration  technique.  The  Gear  integrator  consists  of  a  subroutine 
DIFSUB  which  calls  the  routines  PEDERV,  SOLVE,  and  DIF  FUN.  The  last  routine 

DIF  FUN  must  be  supplied  by  the  user  and  computes  the  rate  of  change  of  the  state 

vector  as  a  function  of  time  and  the  current  values  of  the  state  vector.  For  the 
case  of  equation  (1)  the  state  vector  consists  of  the  variables  z  and  z.  The  rate  of 
change  of  z  is  z  and  the  rate  of  charge  of  z  is 

*  =  m  (M“  hi  ~  **)  (2) 

Equation  (l)  can  be  solved  analytically  if  the  function  f[t)  is  sinusoidal.  Let  us 
consider  the  case 

/(<)  =  (3) 

We  can  assume  a  solution  of  the  form 
•Contributed  by  David  A.  Arnold,  SAO 


which  gives  the  characteristic  equation 


-  AuPm  +  Aiu/b  +  Ak  =  F 


(5) 


Solving  equation  (5)  gives 


A 


F 

-  mu/2  +  iu/b  +  k 


(6) 


If  there  is  no  damping  so  that  b  =  0,  we  have  the  amplitude  given  by 


A 


F 

—  mu/2  +  k 


(7) 


The  amplitude  goes  to  infinity  for 


(8) 


where  Wo  is  the  resonant  frequency.  For  very  high  frequencies  k  can  be  neglected 
and  we  have  a  free  oscillation  given  by 


A  = 


F 

mu/2 


(9) 


For  very  low  frequencies,  mu2  can  be  neglected  giving 
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A 


,  (10) 


At  the  resonance  frequency  given  by  equation  (8),  equation  (6)  becomes 


A 


_F_ 

tub 


(11) 


In  this  case  the  response  is  90°  out  of  phase  with  the  driving  force.  In  equation 
(9)  the  response  is  180°  out  of  phase,  and  in  equation  (10)  the  response  is  in  phase. 
In  addition  to  the  forced  oscillation  given  by  equation  (4)  we  can  have  an 
oscillation  at  the  natural  frequency  given  by 

i  =  Be ■**  (12) 

In  the  case  of  small  damping,  is  given  by  equation  (6).  In  the  presence 
of  damping  the  motion  given  by  equation  (12)  will  decay  with  time  if  there  is  no 
continuing  excitation  of  this  frequency. 

In  the  case  of  a  few  sources  on  a  rotating  wheel  the  function  f(t)  is  not  a 
perfect  sinusoid.  As  the  number  of  sources  is  increased,  the  driving  force  becomes 
more  sinusoidal  but  the  amplitude  becomes  smaller  since  the  distribution  around 
the  wheel  becomes  more  uniform. 

Various  test  runs  have  been  done  with  program  DETECT  to  verify  that  the 
computer  code  gives  the  correct  answer  to  cases  that  can  be  solved  analytically. 
In  the  initial  version,  the  values  of  m,  b ,  k  and  Zb  are  input  parameters.  The 
damping  b  is  specified  as  a  fraction  of  the  critical  damping  coefficient  given  by 
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k  =  2  y/ric  (13) 

A  test  run  with  m  =  k  =  Xe  =  1  and  b  =  0  gave  a  sinusoidal  oscillation  with 
period  2r.  With  m  =  k  =  to  =  1  and  b  =  b^  the  motion  in  a  decaying 
exponential.  The  next  test  run  was  done  with  m  =  k  =  1,  b  =  .01  be,  ip  =  0, 
and 

/(<)  =  /sin  (14) 

with  F  =  1.,  and  =  1  as  computed  from  equation  (8).  The  response  as  a 
function  of  time  is  shown  in  Figure  1.  The  minimums  of  the  curve  occur  at  t  = 
2r,  4  7T,  6tt,  etc.  The  zero  crossings  occur  a  quarter  of  a  cycle  later  when  equation 
(14)  has  its  maximum  value.  This  is  the  phasing  required  to  pump  the  oscillation 
in  the  most  effective  manner. 

The  next  run  was  done  with  m  =  1,  6  =  k  =  Zc  =  0,  and  /  =  1  in 
equation  (14).  The  z  coordinate  increased  monotonically  with  the  velocity 
oscillating  between  0  and  2.  In  equation  (9)  the  amplitude  is  180°  out  of  phase 
with  the  driving  force.  The  value  of  io  has  been  added  to  the  input  parameters 
and  the  run  repeated  with  io  =  -  1.  This  produced  the  results  given  by  equation 
(9)  with  A  =  -  1. 

The  next  test  run  has  been  done  with  m  =  i  =  /  =  l,zc  =  zo  =  0,  and 
b  =  be  =  2.  Using  these  values  in  equation  (6)  with  u>  =  1  gives 

_1_ 

2i 


A 


(15) 
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Figure  2  shows  the  response  as  a  function  of  time.  The  minimums  of  the  curve 
are  at  2t,  At,  6n  etc.  and  the  amplitude  after  the  initial  transient  is  .5  in 
agreement  with  equation  (15).  The  phasing  is  the  same  as  in  Figure  1.  Pumping 
the  oscillator  at  the  resonant  frequency  results  in  a  buildup  of  the  response  until 
the  damping  forces  equal  the  pumping  forces.  From  equation  (4)  we  have  the 
velocity  given  by 


x  =  Aiuc**  (16) 

With  A  given  by  (15)  and  u  -  1,  we  have  i  =  1/2.  Since  b  =  2,  the  damping 
force  is  bx  =  2  x  .5  =  1,  which  is  equal  to  the  driving  force  F  =  1. 

In  the  next  run,  the  damping  b  is  set  to  sero  and  the  other  parameters  are 
kept  the  same  (m  =  k  =  F  =  u  =  l,  and  Xo  -  ip  =  0).  The  results  are  shown 
in  Figure  3.  The  amplitude  is  a  little  larger  than  in  Figure  1  which  had  6  =  .01 

it. 

All  the  runs  so  far  have  had  u>  =  ub.  If  u  is  not  equal  to  u*>  then 
it  is  possible  to  have  components  at  both  frequencies.  The  parameter  u>  has  been 
added  to  the  input  and  set  equal  to  1.5  for  the  next  run.  The  other  parameters 
are  the  same  and  the  natural  frequency  is  m>  =  1.  The  results  are  shown  in 
Figure  4.  The  beats  between  u  and  u*,  can  be  seen.  Since  b  is  sero,  the 
oscillation  at  frequency  ut,  is  undamped. 

At  this  point  some  tests  were  done  to  see  if  the  integration  speed  could  be 
increased  by  decreasing  the  accuracy  requested  from  the  numerical  integrator. 
Reducing  the  accuracy  from  13  to  6  decimal  phases  increased  the  stepsize  from 
.036  seconds  to  .28  seconds.  With  9  decimal  places  the  stepsize  was  .1  seconds. 


Page  9 


With  the  test  runs  completed  using  a  sinusoidal  driving  force  for  the 
detector,  the  program  has  been  modified  to  use  the  signal  input  from  sources  on  a 
rotating  wheel.  The  mathematical  part  of  program  FLYWHEEL  has  been  put  into 
a  subroutine  called  WHEEL  which  has  been  added  to  program  DETECT.  An  input 
variable  specifies  the  number  of  sources  to  be  placed  symmetrically  around  the 
wheel.  As  a  first  test  of  the  program,  a  run  has  been  done  with  m  =  .0001,  b  = 
.1  be,  k  -  u  =  1,  it,  =  ic  =  0,  one  source,  and  the  constant  F  in  the  equation 

FLUX  =  F/r*  (17) 

set  equal  to  625.  The  sources  are  rotating  at  a  45  cm  radius  with  the  detector  at 
70  cm  from  the  center  of  the  wheel.  This  gives  25  cm  for  the  closest  approach 
and  115  cm  for  the  furthest  distance.  From  equation  (17)  the  corresponding  fluxes 
are  1.00000  and  .04726.  The  natural  frequency  i*  from  equation  (8)  with  1=1 
and  m  =  .0001  is  100  hertz.  This  gives  a  fast  rise  time  for  the  detector  so  that  it 
will  record  the  instantaneous  value  of  the  flux  for  test  purposes.  Figure  5  shows 
the  response  of  the  detector  to  the  rotating  source.  The  maximum  response  is 
1.00  and  the  minimum  is  .04726  in  agreement  with  the  calculated  values  of  the 
flux.  The  curve  shape  looks  the  same  as  curve  shapes  generated  by  program 
FLYWHEEL 

In  the  actual  experiment  the  rotating  wheel  will  spin  so  as  to  produce  a 
signal  having  the  same  frequency  as  the  <4,  of  the  detector.  The  next  run  has 
been  done  with  1  source,  and  k  =  m  =  1  to  produce  a  natural  frequency  14,  =  1. 
The  angular  velocity  of  the  wheel  is  u  =  1,  and  carries  one  source.  The  d*mpin£ 
b  is  set  to  zero.  Figure  6  shows  the  response  of  the  detector.  The  input  signal  is 
always  positive  so  that  the  mean  position  has  a  slight  positive  base. 
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In  order  to  compare  the  response  of  the  detector  to  i  otating  source  with 
the  response  to  a  sinusoid  of  the  same  amplitude,  a  run  has  been  done  with 
a  sinusoid  of  amplitude  F  =  .5.  The  peak  to  peak  driving  force  is  1  with  the 
signal  varying  between  -.5  and  4.5.  For  the  case  of  a  single  rotating  source 
the  peak  to  peak  amplitude  is  .95274  and  the  signal  varies  between  .04726  and 
1.00000.  Figure  7  shows  the  response  to  a  sinusoid  of  amplitude  F  =  J>. 

A  series  of  runs  have  been  done  with  2,  3,  4,  and  5  sources  uniformly 
spaced  around  the  wheel.  The  rotation  frequency  of  the  wheel  is  1/n  where  n  is 
the  number  of  sources.  The  response  of  the  detector  is  shown  in  Figures  8,  9,  10 
and  11  respectively. 

The  phasing  of  the  rotation  of  the  wheel  is  such  that  the  first  source  is  at 
its  closest  approach  to  the  detector  at  t  =  0,  thereby  giving  a  maximum  in  the 
signal  at  t  =  0.  This  is  90°  out  of  phase  with  the  signal  given  by  equation  (14). 
For  the  sake  of  easier  comparison  of  results  the  program  has  been  changed  to  use 
the  equation 


f(t)  =  F  sin  (u*>t  +  6)  (18) 

where  6  is  a  phase  angle.  Setting  6  =  90°  causes  f{t)  to  have  a  maximum 
at  1  =  0  the  same  as  the  signal  from  the  rotating  wheel.  A  run  has  been  done  to 
give  the  response  of  the  detector  using  F  *  .5,  6  =  90°,  m  *  k  =  u  =  1,  and 
2c  =  io  =  0.  The  result  is  shown  in  Figure  12. 

Table  1  gives  a  summary  of  the  response  of  the  detector  to  various  types 
of  signals  as  shown  in  Figures  6  through  12. 
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Table  1 


Figure 

Signal 

Time 

*(0 

Ax 

6 

1  source 

45.586 

48.728 

♦  5.895 
-5.792 

11.687 

7 

.5sinu/t 

40.840 

43.982 

47.124 

4  10.210 

•10.996 

411.781 

21.206 

22.777 

8 

2  sources 

45.609 

48.750 

48.390 

-8.058 

16.448 

9 

3  sources 

45.635 

48.774 

48.975 

-8.315 

17.290 

10 

4  sources 

45.668 

48.805 

48.541 

-7.480 

16.021 

11 

5  sources 

45.715 

48.848 

7.670 

-6.166 

13.836 

12' 

.5cosurt 

45.574 

48.715 

4  11.391 
-12.176 

23.567 

Table 

1.  Response  of  the  detector 

to  various 

input  signals. 

For  all  of  the  cases  except  Figure  7,  there  is  a  minimum  at  about  45.6  seconds 
and  a  maximum  at  around  48.8  seconds.  The  last  column  lists  the  difference 
between  the  peaks.  Figures  7  and  12  which  are  sinusoidal  signals  give  the  largest 
response.  For  the  case  of  n  sources  on  a  rotating  wheel,  the  signal  is  more 
sinusoidal  as  the  number  of  sources  is  increased,  but  the  difference  between  the 
peaks  of  the  driving  signal  decreases  as  the  pulses  from  individual  sources  start  to 
overlap.  The  largest  output  response  in  Table  1  occurs  for  3  rotating  sources. 

As  a  measure  of  the  efficiency  of  the  various  types  of  driving  signals 
the  ratio  of  the  detector  response  to  the  signal  amplitude  has  been  computed 
for  each  of  the  signal  types  listed  in  Table  1  (except  for  Figure  7  which  is 
difficult  to  compare  because  the  phase  is  different).  Table  2  lists  the  type  of 
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source,  the  minimum  value  of  the  signal  Fm in.  the  maximum  value  the  peak 
to  peak  amplitude  AF,  the  peak  to  peak  output  response  Ax  taken  from  Table  1, 
the  ratio  of  the  output  to  the  input  Ax/ AF,  and  the  efficiency  relative  to  .5  cos  ut 
computed  by  dividing  the  values  of  Ax/AF  by  the  last  entry  in  that  column. 

Table  2 


Signal 

F twin 

F mxx 

AF 

Ax 

Ax 

AF 

Efficiency 

1  source 

.047 

1.000 

.953 

11.687 

12.266 

.52 

2  sources 

.152 

1.047 

.895 

16.448 

18.368 

.78 

3  sources 

.303 

1.114 

.811 

17.290 

21.317 

.90 

4  sources 

.494 

1.199 

.705 

16.021 

22.708 

.96 

5  sources 

.710 

1.300 

.590 

13.836 

23.428 

.994 

.5cosut 

-.500 

.500 

1.000 

23.567 

23.567 

1.000 

Table  2. 

Efficiency  of 

various 

input  signals  in 

driving  the 

detector 

The  detector  response  has  been  calculated  after  a  fixed  interval  of  time, 
since  there  is  a  resonant  buildup  in  the  amplitude.  Since  no  damping  was  used, 
the  buildup  would  continue  indefinitely.  Equation  (6)  gives  the  amplitude  for  any 
frequency  and  with  damping  in  the  system.  Some  runs  have  been  done  to  test  the 
program  under  other  conditions.  For  example,  for  m  =  k  -  u*,  =  F  =  1,  b  =  0, 
and  w  =  1.5  equation  (6)  gives  A  =  -  .8.  The  detector  is  being  driven  above  its 
resonant  frequency,  and  the  response  is  180°  out  of  phase  with  the  driving  force. 
A  run  has  been  done  with  Xo  =  -  .8  and  ic  =  0.  The  results  are  shown  in  Figure 


Tirnp 


-li. 


Timr 


Time  (see) 


Time 
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13.  The  response  is  a  steady  state  oscillation  as  described  by  equation  (6). 

If  the  damping  is  b  =  .lit  with  m  =  k  =  1,  the  value  of  b  from  equation 
(13)  is  .2.  If  u>  =  wo  =  F  =  1,  the  response  from  equation  (11)  is  A  =  5/t.  A 
run  has  been  done  with  Xo  =  ie  =  0.  The  results  are  shown  in  Figure  14.  After 
an  initial  transient  phase,  the  amplitude  approaches  8  asymmetrically.  The  last 
peak  has  an  amplitude  of  -4.96.  The  driving  force  is  a  cosine  function  which  has 
zero  phase  at  6.283,  12.566,  18.849,  25.133,  31.416,  37.699,  and  43.982  seconds. 
Toward  the  end  of  the  run,  Figure  14  has  maxima  at  32.984,  39.267,  and  45.551 
seconds,  which  is  jt/2  seconds  after  the  maxima  of  the  driving  function. 

Finally  a  run  has  been  done  with  u>  =  .5,  m  =  k  =  Uo  =  F  =  l,  and 

5  =  0.  From  equation  (6)  we  have  A  =  4/3.  The  detector  is  being  driven 

« 

below  its  resonant  frequency  ^nd  the  response  is  in  phase  with  the  driving  force. 
A  run  has  been  done  with  Zo  =  4/3,  and  io  =  0.  The  response  is  shown  in 
Figure  15.  The  behavior  is  exactly  as  given  by  equation  (6). 


3.0  SIGNAL  FROM  3  ROTATING  SOURCES  VS.  WHEEL  DIAMETER 

From  an  engineering  point  of  view,  it  is  desirable  to  keep  the  rotating  wheel 
carrying  the  neutrino  sources  smaller  in  order  to  minimize  centrifugal  effects.  A 
series  of  runs  have  been  done  for  different  radii  to  see  the  effect  on  the  signal  from 
three  evenly  spaced  sources.  It  is  assumed  the  closest  approach  to  the  detector  is 
25  centimeters  for  all  the  cases.  The  flux  is  computed  using  equation  (17)  with 
F  =  625.  This  gives  unity  at  25  cm.  Table  3  gives  the  results  of  the 
calculations.  The  first  column  is  the  radius  of  the  wheel,  the  second  is  the 
maximum  signal,  the  third  is  the  minimum  signal,  and  the  last  column  is  the  peak 
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to  peak  amplitude  of  the  signal.  As  the  wheel  gets  smaller,  the  peak  to  peak 
amplitude  is  reduced  but  the  dependence  is  not  severe  over  the  range  of  values 
studied. 


Table  3 


r(em) 

Ftn'm 

AF 

34.5 

1.1717 

.4519 

.7199 

40.0 

1.1374 

.3638 

.7736 

43.0 

1.1228 

.3256 

.7972 

45.0 

1.1143 

.3033 

.8111 

51.5 

1.0921 

.2438 

.8482 

Table  3.  Neutrino  signal  amplitude  vs.  wheel  radius  using  three  symmetrically 
placed  sources. 


4.0  GRAVITY  SIGNAL  VS.  RADIUS  OF  THE  ROTATING  WHEEL 

The  previous  section  discusses  the  effect  of  wheel  radius  on  the  neutrino 
signal  from  three  rotating  sources.  In  order  to  minimize  the  gravity  signals  from 
the  wheel,  dummy  sources  consisting  of  empty  9  cm  spheres  are  placed  between 
the  neutrino  sources.  As  the  wheel  radius  is  decreased,  fewer  balls  will  fit  on  the 
rim  and  the  gravity  fluctuation  increases  as  the  wheel  radius  decreases.  The 
gravity  flux  has  been  computed  using  equation  (17)  with  a  coefficient  of  625  so 
that  the  signal  is  unity  for  one  source  at  25  cm.  Table  4  gives  the  results  of  the 
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computations.  Tbe  first  column  is  the  wheel  radius,  the  second  is  the  number  of 
masses,  the  third  is  the  maximum  gravity  signal,  the  fourth  column  is  the 
minimum,  and  the  fifth  is  the  peak  to  peak  amplitude.  The  absolute  magnitude 
of  the  signal  is  listed  (the  actual  force  being  negative).  With  a  45  cm  radius,  it  is 
not  possible  to  fit  36  nine  centimeter  diameter  balls  since  the  circumference  is 
282.7  cm  and  324  cm  is  needed  for  36  balls.  For  all  the  other  cases,  the  balls 
should  be  just  touching.  For  the  45  cm  radius  the  balls  would  have  to  be  7.85  cm 
instead  of  9  cm.  This  slightly  closer  spacing  gives  a  significant  reduction  in  the 
gravity  ripple.  However,  the  gravity  ripple  is  a  factor  of  104  lower  than  the  signal 
from  a  single  ball  for  the  first  entry  which  is  the  largest  one  in  the  table.  The 
magnitude  of  the  signal  from  all  the  balls  is  equivalent  to  about  6  balls  at  the 
closest  approach  (25  cm  from  the  detector). 

Table  4 


r(cm ) 

#masses 

fni 

nun 

AF 

34.5 

24 

5.8068743 

5.8067482 

.0001261 

40 

28 

5.9817155 

5.9816336 

.0000819 

43 

30 

6.0197642 

6.0196919 

.0000723 

45 

36 

6.9407545 

6.9407441 

.0000104 

51.5 

36 

6.1496537 

6.1496060 

.0000477 

Table  4.  Gravity  signal  vs.  wheel  radius  and  number  of  masses. 

For  the  first  entry,  the  signal  is  produced  by  24  balls,  three  of  which  are 
neutrino  sources.  The  frequency  of  the  gravity  signal  is  8  times  that  of  the 
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neutrino  signal.  Let  us  assume  the  detector  is  tuned  to  the  frequency  of  the 
neutrino  sourtes.  If  we  neglect  damping,  we  can  use  equation  (7)  to  calculate  the 
off-resonance  response  of  the  detector.  For  k  =  m  =  F  =  1,  and  u  =  8, 
equation  (7)  gives  A  =  .01587.  For  u  =  0,  the  amplitude  from  either  equation  (6) 
or  (7)  is  1.  At  resonance  the  amplitude  is  given  by  equation  (11).  If  we  have  a 
high  Q  system  the  damping  coefficient  b  given  by  b  -  mu/Q  will  be  small  and 
the  amplitude  A  will  be  large.  For  u  =  1,  and  Q  =  104,  we  have  b  =  10~4.  The 
amplitude  from  equation  (11)  is  A  -  -  104!.  The  amplitude  at  u  =  8  is  about  2 
orders  of  magnitude  smaller  than  the  response  at  u;  =  0,  and  about  6  orders  of 
magnitude  smaller  than  the  response  at  resonance  with  Q  =  104. 

Up  to  this  point,  the  detector  has  been  modelled  as  a  one  degree  of  freedom 
harmonic  oscillator.  This  is  probably  a  reasonable  approximation  for  the  neutrino 
signal  on  the  sapphire  crystal.  However,  for  the  effect  of  gravity,  the  acceleration 
of  the  frame  must  be  considered.  In  addition  to  translation,  the  frame  can  also 
rotate  due  to  unequal  force  on  the  upper  and  lower  sensing  masses. 


5.0  TWO  MASS  MODEL  OF  THE  NEUTRINO  DETECTOR 

The  neutrino  detector  consists  of  a  sensor  plate  that  is  attached  to  a  frame. 
The  frame  is  suspended  by  a  soft  suspension  system  to  minimize  the  transmission 
of  seismic  noise  to  the  detector.  An  analysis  has  been  done  of  the  response  of  the 
system  using  the  model  shown  in  Figure  16.  The  equations  of  motion  of  the 
system  are 
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Figure  16.  Two  mass  model  of  the  neutrino  detector  with  damping. 


miij  +  b\i\  -  hfa  -  xi)+  ( ki  +  *:)ii  -  kjZi  =  Fie*wt  (19) 

mjfj  +  fcj(ij  -  *i)+  kin  -  ktzi  =  Fte**  (20) 

Assuming  the  solutions 

*l  =  Axt **  (21) 

Z7  =  A,e~‘  (22) 


gives  the  equations 

-  Aimju;1  +  A\ib\u>  -  »7>jw(Aj  -  Aj)  +  [k\  +  ki)A\  -  k^Ai  ~  F\  (23) 

-  Ajtoju;*  +  ib2u{Ai  -  Ai)+  AjAj  -  ijAl  =  Ft  (24) 

Equations  (23)  and  (24)  can  be  written  in  the  form 


A\B  +  AjC  =  F\ 
A\C  +  AtD  =  Fi 


(25) 

(26) 
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where 


B  =  -  miJ1  +  +  kt  +  iu;(5i  +  bt) 

(27) 

C  =  -  kt  -  ibtu 

(28) 

D  =  -  mtu?  +  kt  +  ibtu> 

(29) 

Equations  (25)  and  (26)  can  be  solved  to  give 


Ai  =  Xi/A  (30) 

At  =  Xt/A  (31) 

where 

Xi  =  FiD-FtC  (32) 

X,  =  FtB-FiC  (33) 

A  =  BD-C*  (34) 


In  order  to  evaluate  equations  (32),  (33),  and  (34)  let  us  write  equations  (27),  (28) 
and  (29)  in  the  form 


B  =  Bi  +  »  Bt  (35) 

C  *  Ci  +  i  Ct  (36) 

D  ■  D\  +  «  Dt  (37) 

where 

B\  =  +  +  (38) 

Bt  =  u{bi  +  h)  (39) 

Ci  m  -kt  (40) 
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Ci  = 

-  b2uj  =  -  Di 

(41) 

Di  = 

-  mju>*  +  k2 

(42) 

Di  = 

hi  U) 

(43) 

Substituting  equations  (35)  through  (37)  into  equations  (32)  and  (33)  gives 

Xi  =  F\D\  -  Ft Ci  +  i(FiDi  -  FiCh)  (44) 

Xt  =  FtBi  -  FiCi  +  .•  [FtBt  -  FiCt)  (45) 

Substituting  equations  (38)  through  (43)  into  (44)  and  (45)  gives 

X\  =  mjws  +  kij  +  /jij  +  «(fi  +  (46) 

Aj  =  /j(- miw*  +  ii  +  + +  »w^jF:[5i  +  5j]  + (47) 

Substituting  equations  (35)  through  (37)  into  equation  (34)  gives 

A  =  (Bi  +  iBi)  (l>i  +  «£>:)  -  (C,  +  »C2)* 

=  BiDi  -  B2Di  -  Ci*  +  C:*  +  »(b:Z?i  +  £i/?2  -  2CiCj)  (48) 


We  can  write  equation  (48)  as 


^  =  As  +  I 


(49) 


where 


4a  =  -  B2Di  -  Ci1  +  Cf 

A i  =  BjDi  +  B\Di  -  2C\Ci 


(50) 

(51) 
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Evaluating  the  terms  in  equation  (50)  gives 


(52) 

(53) 

(54) 

(55) 

Substituting  equations  (52)  through  (55)  into  (50)  gives 


BiDi  =  mjmjw4  -  w*  (miij  +  mjiti  +  mjij)  +  *i*j  +  *2: 

-  BtDi  —  - 

=  -  *1 
Cj1  = 


4*  =  mimju;4  -  miiks  +  mjiki  +  mjij  +  5i6j^  +  tiikj  (56) 

Evaluating  the  terms  in  equation  (51)  gives 

BjDi  =  -  4- 

BiDj  =  —  +  hiki'j 

2C]C2  s  —  2u>&2&2 

Substituting  equations  (57)  through  (59)  into  (51)  gives 

4/  =  -  w3(mi5j  +  ro2&i  +  +  u/^Mj  +  Ml)  (60) 

Combining  equations  (49),  (56)  and  (60)  gives 


(57) 

(58) 

(59) 


4  =  numju4  -  vi^miki  +  miki  +  miki  +  bih'j  +  k\kt 
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+  “  wS(fni&*  +  +  mjfcjj  4  u(^bik*  4  bik\  jj  (61) 

Equations  (46)  and  (47)  can  be  written  in  the  form 


X\  =  Xir  +  i  Xu  (62) 

Xt  =  Xir  4  i  Xu  (63) 

where 

Xir  =  Fi(-  mtu/3  +  k:)  +  fj/ki  (64) 

Xu  =  (ft+Ftjhv  (65) 

Xt*  =  -  miu>*  +  kt  +  kt j  +  Fiki  (66) 

Xu  =  4  6j]  +  Fib^j  (67) 


Substituting  equations  (49),  (62),  and  (63)  into  equation  (30)  and  (31)  gives 


j4i  =  +  iXu 

Ar  4  iAj 

Aj  =  ±  jXu 

Ar  4  iAj 


(68) 

(69) 


Equations  (68)  and  (69)  can  be  rationalised  by  multiplying  numerator  and 
denominator  by 


A*  =  Ar  —  iAi 


(70) 


The  denominator  becomes 
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d*  —  dd*  —  dj?  +  A / 

pi) 

Multiplying  equations  (62)  and  (63)  by  (70)  gives 

XxA'  =  (Xut+iXu)  (d*  -  id/) 

(K) 

XtA *  =  ( XtM+iXti )  (ar  -  id/) 

PS) 

or 

•Kid*  =  (A'ued*  Jfi/d/)  +  i^A’i/d/e  -  A'ujd/) 

(74) 

XiA'  =  ^i/jdi  +  A’j/d/)  +  i^Xt/A/i  -  XtpAjj 

(75) 

The 

final  results  are 

A/  *  XiA‘/A* 

(76) 

Aj  =  XjA'/A2 

(”) 

with 

the  terms  defined  by  equations  (71),  (74)  and  (77). 

If  Ar  is  zero,  equations  (68)  and  (69)  simplify  to 

^  _  ATxi?  +  iXu 
iAj 

(78) 

l  _  XtR  +  iXu 
id/ 

(79) 

or 

Ai  =  Xu/Aj  -  iXir/Aj 

(80) 

M  =  Xu/Aj  -  iXt r/Aj 

(81) 
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In  the  special  case  of  6]  =  0,  we  see  from  equations  (65)  and  (67)  that  Xu  =  Xu 
so  that  A\  and  Aj  have  the  same  real  component  and  differ  only  in  the  imaginary 
part.  In  general  the  amplitudes  are  given  by  equations  (76)  and  (77). 

In  the  case  of  a  high  Q  system,  Aj  will  have  a  small  value.  The  amplitude 
of  the  response  will  be  very  large  when  Ar  goes  to  zero.  In  the  case  of  zero 
damping  we  see  from  equation  (60)  that  Ai  is  zero,  and  the  amplitudes  Ai  and  A2 
will  go  to  infinity  when  Ar  is  also  zero.  Setting  equation  (56)  to  zero  gives  two 
frequencies  u>i  and  wj  where  Ar  =  0.  These  are  the  resonant  frequencies.  For 
small  damping  (high  Q)  the  maximums  of  Ai  and  A}  are  nearly  at  wi  and  wj  so 
equations  (80)  and  (81)  can  be  used  to  compute  the  maximum  amplitudes. 

Some  case  studies  have  been  using  the  parameters 

mi  =  12.5  kg 
mj  =  .25  kg 

k\  ss  500  newtons /meter  (82) 

kt  =  10s  newtons / meter 

For  systems  1  and  2  individually  the  resonant  frequencies  given  by 


are  ui  —  6.324555320  rad  /sec  and  =  632.455320  rad/sec.  Using  a  value  of 
h  =  .016  newtons/(m/sec)  in  the  equation 


Q  =  mu/b 


(84) 
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gives  a  value  of  9862  for  Qi,  or  about  104.  Using  a  value  of  bi  =  4 
newtons/(m/sec)  gives  a  value  of  19.76  for  Q\,  or  about  20.  Using  the 
parameters  given  by  equation  (82)  with  no  damping  the  roots  of  equation  (56)  are 
u>i  =  6.262236891  rad/sec  and  wi  =  638.7493909  rad/sec.  The  first  frequency 
corresponds  to  masses  mi  and  mi  moving  together.  The  frequency  is  somewhat 
lower  than  for  system  1  alone  because  the  mass  is  larger.  The  second  frequency 
corresponds  to  masses  mi  and  mi  moving  opposite  to  each  other.  The  frequency 
is  somewhat  higher  than  for  system  2  alone  because  the  oscillation  is  about  the 
center  of  mass  and  the  effective  spring  constant  is  higher  than  k i. 

Four  cases  have  been  run  with  the  parameters  given  by  equation  (82)  plus 


1)  s  0 

bi  =  .016  (85) 


for  the  first  3  cases  and  bi  =  4  for  the  fourth. 


The  cases  differ  in  the  type  of  applied  forces  F\  and  Ft.  All  the  cases  are 
at  the  second  resonance  given  by  u>j. 


Case 


1 

2 

3 


Fi 

Ft 

0 

(bi=0) 

m2 

mi 

.86mi 

mi 

mi 

0 

(b,=4) 

mi 

Amplitude 


Ai  =  -  .4806  x  10-7  +  .4703  x  lO"3* 
Ai  =  -  .4806  x  10-7  -  .2351  x  10"1* 

Ai  =  -  .2444  x  10-*  -  .6589  x  10"4i 
At  =  -  .2444  x  10-*  +  .3294  x  10-*i 

A i  =  -  .2451  x  10-*  -  .4611  x  10“7i 
A :  =  -  .2451  x  10-*  +  .2305  x  10"5* 

Ai  =  -  .4385  x  10-7  +  .4290  x  lO"3* 
Ai  =  -  .1101  x  10-4  -  .2145  x  10-1* 


4 
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For  cases  1,  2,  and  3,  the  real  parts  of  Ai  and  A:  are  equal  (see  discussion 
following  equations  (80)  and  (81)).  For  all  except  case  3,  the  real  part  of  the 
amplitude  is  small  compared  to  the  imaginary  part.  The  amplitudes  are  limited 
by  the  damping  and  the  response  is  90°  out  of  phase  with  the  driving  force.  For 
case  3  the  real  and  imaginary  parts  of  At  are  nearly  equal  and  the  phase  angle  is 
134  degrees.  For  all  4  cases  the  ratio  of  the  imaginary  part  of  At  to  the 
imaginary  part  of  A\  is  *50.  The  ratio  of  the  masses  is  12.5/.25  =  50.  The 
motions  of  mi  and  mj  are  inversely  proportional  to  their  masses  and  180°  out  of 
phase  (except  for  case  3).  If  we  normalize  the  imaginary  parts  of  At  to  that  of 
case  1,  the  relative  amplitudes  of  the  imaginary  parts  are 

Case  Amplitude 

1  l.Od 

2  .14 

3  .0001 

4  .91 

The  relative  amplitudes  of  the  4  cases  are  determined  primarily  by  the  difference 
in  the  driving  accelerations  between  the  masses.  In  the  first  case,  all  the 
acceleration  is  applied  to  the  second  mass.  This  case  gives  the  largest  response  of 
the  system.  In  the  second  case  the  acceleration  difference  is  14  percent  of  the 
acceleration  difference  in  case  1.  The  response  is  lower  by  14  percent.  In  case  3 
the  applied  accelerations  are  equal.  This  gives  the  smallest  response.  The 
difference  A\  -  At  is  not  sero  because  of  the  effect  of  k\.  For  a  single  mass, 
the  equation  of  motion  is 


mi  •+  bz  +  k  =  Ft"* 


(86) 


Assuming  a  solution  of  the  form 
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x  =  Ac* 


results  in  the  solution 


—  thu*  +  t  +  iuv 

At  resonance,  the  real  part  of  the  denominator  goes  to  zero  and  the  response 
is  90°  out  of  phase  with  the  driving  force  F.  Above  the  resonant  frequency, 
the  term  -  m  J1  is  much  larger  than  k.  For  the  4  cases  that  have  been  studied 
the  driving  force  is  close  to  the  natural  frequency  of  mj  but  is  about  100  times  the 
resonant  frequency  of  mi  on  spring  k\.  As  free  masses  {b  =  k  =  0)  with  an 
applied  force  ft,m,  the  response  of  each  mass  using  equation  (88)  is 

A  ■  ■‘fH  ~  *  (89> 

-  mw*  ur 

If  *2  is  present  but  k\  is  zero,  equation  (89)  still  applies  since  Aj  -  Ai  is  always 
zero.  Finally,  adding  spring  k\  causes  a  small  force  on  A\  which  introduces  a 
small  difference  between  A\  and  A}.  The  approximate  magnitude  of  the 
perturbation  caused  by  k\  can  be  calculated  as  follows.  At  the  resonant  frequency 
u>i  of  mass  1,  we  have 

-  miwi*  =  k\  (90) 


At  frequency  where  wj/u>i  =  105,  we  have 
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to  the  line  width  of  about  .1  rad /sec  at  half  amplitude. 

Figures  17,  16,  and  19  show  the  resonance  plots  for  cases  1,  2,  and  3 
respectively.  The  plots  for  case  4  look  very  similar  to  case  1  in  Figure  17  and 
have  been  omitted.  The  quantities  plotted  in  parts  a),  b),  and  c)  are  the  r*?l  part 
of  Aj,  the  imaginary  part  of  Aj  and  the  magnitude  of  Aj.  Because  of  the  unusual 
behavior  of  Figure  19c,  the  magnitude  has  bee*,  plotted  over  a  wider  region  in 
Figure  19d.  The  last  part  of  each  figure  is  the  phase  angle.  Figures  17  and  18 
for  cases  1  and  2  look  like  standard  resonance  plots.  For  case  1  the  real  part  is 
positive  below  resonance  and  negative  above  resonance.  The  imaginary  part  and 
the  magnitude  have  their  maximum  values  at  resonance.  The  phase  angle 
approaches  0°  below  resonance  and  approaches  180°  above  resonance.  For  case  2, 
the  difference  in  applied  force  is  of  opposite  sign  and  the  displacements  are 
reversed  in  sign  and  phase.  Case  3  is  unlike  all  the  others.  For  the  real  part  is 
always  negative.  The  shape  of  the  curve  is  similar  to  that  of  case  2.  The 
imaginary  part  is  also  similar  to  case  2.  The  major  difference  is  in  the  magnitude. 
There  is  no  resonance  peak.  The  shape  is  similar  to  the  real  part  of  the  amplitude 
in  case  1.  The  maximum  amplitude  occurs  below  resonance.  The  value  at 
resonance  is  comparable  to  that  on  either  side  of  the  resonance.  The  phase  angle 
approaches  180°  on  either  side  and  decreases  by  about  46°  to  around  134° 
resonance. 

In  summary  the  study  of  the  two  mass  model  shows  that  the  response  at 
the  un  resonance  is  proportional  to  the  difference  in  acceleration  applied  to  the  two 
masses.  The  response  for  equal  accelerations  is  very  small  but  does  not  go  to  zero 
because  of  the  effect  of  the  spring  Aj. 


Figure  19(c). 
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6.0  NUMERICAL  INTEGRATION  OF  THE  TWO  MASS  MODEL 

Section  2.0  describes  a  program  called  DETECT  which  integrates  the  motion 
of  a  single  mass  acted  on  by  an  elastic  restoring  force,  a  damping  force  and  a 
driving  force  which  may  be  sinusoidal  or  the  signal  from  a  rotating  wheel  carrying 
various  numbers  of  rotating  sources.  Section  5.0  gives  the  equations  of  motion  for 
a  two  mass  system  where  one  mass  can  represent  the  frame  and  the  other  mass 
can  represent  the  detector  attached  to  the  frame.  A  ne*-  version  of  program 
DETECT  called  DETECT!  has  been  written  using  the  equations  of  motion  in 
Section  5.0.  The  new  program  can  handle  either  1  or  2  masses.  For  one  mass 
the  program  is  identical  to  program  DETECT.  The  equation  of  motion  for  mass  1 
is 

m\i\  =  -  bii\  -  k\x\  +  Fie1"1 

+  ~  km  +  km  (94) 

If  the  single  mass  mode  is  selected  the  terms  involving  5:  and  ki  are  omitted. 
The  equation  of  motion  for  mass  2  is 

mjij  =  -  5j*j  -  km  +  Fjew‘ 

+  fcxi  +  km  (95) 

The  steady  state  solution  of  equations  (94)  and  (95)  is  given  in  Section  5.0 
by  equations  (61)  through  (69).  The  solution  for  each  of  the  masses  is  of  the 
form 


x 


A  «“*" 


(96) 
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where  the  amplitude  A  is  a  complex  number 

A  =  Ar  +  iAj 

The  driving  force  is  of  the  form 

/  si  pc’**  =  F  cosut  +  iF  tin  urt 
Substituting  equation  (97)  into  equation  (96)  gives 

x  =  (a*  +  iA/)  (co6  ut  +  isinu/t 

At  t  =  0,  we  have 

to  »  An  +  iAi 

The  velocity  obtained  by  differentiating  equation  (96)  is 

x  =  luAe^1 

Substituting  equation  (97)  into  equation  (101)  gives 

x  =  »u;(A*+iA/)  (cosu it  +  iainurt 


(97) 


(98) 


(99) 


(100) 


(101) 


(102) 


At  t  =  0,  we  have 
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io  =  -  uAj  +  iuAjt  (103) 

Taking  the  real  parts  of  equations  (98),  (100),  and  (103)  we  have 

/  =  F  cos  ut  (104) 

to  =  Ar  (105) 

io  -  -uAj  (106) 

If  A  is  real,  the  initial  '-elocity  is  aero  and  the  initial  conditions  consist 
of  giving  only  the  initial  position.  If  A  is  pure  imaginary  the  initial  position 
is  zero  and  there  is  only  an  initial  velocity.  If  A  is  pure  imaginary,  it  is  possible 
to  compute  initial  conditions  involving  only  the  position  by  means  of  the  phase 
shift 


a*  -  ut'  +  | 

cos  ut  =  sin  u >t ' 
sin  ut  =  -  co6  u rt  i 


(107) 

(108) 
(lt») 


Substituting  equations  (108)  and  (109)  into  equations  (98),  (99),  and  (102)  gives 


/=  F  Bin  ut '  -  iF  co6  ut ' 
z  =  ^  Ar  +  iAi  j  (sin  ut '  -  i  co6  ut  '  ) 


(110) 

(HI) 


x  =  iu^Ar  +  iAj^j  (sin  ut  '  -  i  cos  ut ' ) 
=  ( iuAr  -  uA/'j  (sin  ut1  -  i cos  ut ' ) 


(112) 


At  t  =  0,  the  real  parts  of  equations  (111)  and  (112)  are 


*o  «  Ai  (113) 

*»  *  uAr  (114) 

The  real  parts  of  equation  (110)  is 

/  =  Fsinurt'  (115) 

Equations  (113)  and  (115)  can  be  used  in  cases  where  Ar  is  zero. 

Various  test  runs  have  been  done  to  test  the  numerical  integration  program 
DETECT2.  In  the  first  test  a  case  run  or  program  DETECT  was  repeated  on 
program  DETECT2  in  the  1  mass  mode  with  identical  results. 

For  the  two  mass  tests  the  parameters  used  are 


mi 

=  12.5 

(116) 

*1 

=  500 

(117) 

m* 

=  .25 

(118) 

** 

=  10s 

(110) 

If  b\  -  6j  =  k\  =  0,  and  Fi  =  -  F%,  the  equations  of  motion  are  quite  simple. 
Equations  (56),  (64),  (66)  from  Section  5.0  become 

Ar  =  mjmjw4  -  u/J  +  mtkz'j 

=  wJ  -  *:[mj  +  m2]j  (120) 
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X\r  =  -  /*(- mju/-  .ij j  +  Fjki 

=  fjmjw5  (121) 

X\ R  =  Fj(-  miw1  +  A*)  -  /jAj 

=  -  Fjmiu/*  (122) 


The  amplitudes  are 


X\r  _  _ mtFi _ 

Ar  -  Aj(mi  +  mj) 

_ =Jj _ 

-  miv*  +  Aj(mi  +  mj)/mj 


(123) 


XiR  _  _ _ -  *»>/» _ 

4*  -  A*(mi  +  mj) 

; _ it _ 

-  mjw*  +  ki[mi  +  mj)/mi 


(124) 


The  equation  of  motion  for  a  single  oscillation  without  damping  given  by  equation 
(7)  in  Section  2.0  is 


A 


F 

mw5  +  A 


(125) 


Equations  (123)  and  (124)  can  be  written 


A\ 


-Ft 

-  mjw*  +  Aj' 


(126) 
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rl 

-  mtv*  +  kt 

(127) 

where 

ki'  =  ki^mi  +  mjj/mj 

(128) 

(129) 

For  the 

parameters  given  fay  equations  *(116),  (118)  and  (119), 

we  have 

=  51  kt  =  51x10s 

'  (130) 

*»'  *  56  *»  = 

(131) 

The  natural  frequency  for  equation  (125)  is 

• 

V  m 

(132) 

The  natural  frequency  for  either  equation  (126)  or  (127)  is 


<•* 


s 


n»i  +  mj 


=  638.7487769  rad/sec 


(133) 


At  M„,  the  amplitude  is  infinite.  Program  TWOMASS  has  been  used  to  plot  the 
response  from  630  to  640  rad/sec.  At  640  rad/sec,  we  have  the  particularly 


simple  response 

Ax  =  SxlO-6  Ft  (,o-) 

Ai  =  -  2.5  x  10*3  Ft  (135) 

A  computer  simulation  has  been  done  with  program  DETECT 2  using  the 
initial  conditions  *j  =  5  x  10-*,  ij  =  0,  xj  =  -  2.5  x  10-3,  xj  =  0,  Fi  =  -  1, 
fj  s  +  1,  b  «  10s,  =  0,  mi  =  12.5,  mj  =  .25  and  w  =  64C  rad /sec. 

The  driving  function  was 

f\  =  -  cos  urt  (136) 

/j  =  cos  ut  (137) 

The  maxima  of  Xi  were  -  5  x  10-s  at  .00490874  sec  and  5  x  10-s  at  .00981747. 
With  a  period  of  .00981748,  the  maxima  are  at  .5  and  1.0  cycles.  The  maxima  of 
x;  were  2.5  x  10_l  at  .00490874  seconds  and  -  2.5  x  10~3  at  .00981747  seconds. 
The  agreement  with  the  theoretical  calculations  was  essentially  exact  within  the 
limits  of  the  interpolation  error  between  output  points. 

Program  TWOMASS  has  been  rerun  with  the  same  conditions  as  above  but 
adding  5}  =  32.  This  gives  about  10%  of  critical  damping.  The  amplitudes 
computed  by  the  program  are 

A\  -  0  +  9.5928966  x  10-7  i  (138) 

At  =  0  -  4.7964483  x  10- 8  »  (139) 

Program  DETECT2  has  been  run  with  the  initial  values  of  xi  and  xj  set  equal  to 
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the  imaginary  parts  of  A\  and  At.  The  initial  velocities  are  zero  and  the  driving 
function  is 


f\  =  -  sin  ut  (140) 

ft  =  sinu/f  (141) 


The  initial  conditions  and  driving  force  are  described  by  equations  (113),  (114),  and 
(115).  The  driving  frequency  is  the  resonant  frequency  given  by  equation  (133). 
The  period  is  .009836708  seconds.  The  maxima  of  x\  were  -9.5928960xl0*7  at 
.00491835  seconds  and  9.5928943 xlO*7  at  .00983670  seconds.  The  maxima  of 
z:  were  4.7964480x10-*  -at  .00491835  seconds  and  -4.796447 lxlO*5  at  .00983670 
seconds.  The  amplitude  at  resonance  for  a  single  oscillation  with  damping  is 


A 


F_ 

ibu 


(142) 


as  shown  in  equation  (11)  of  Section  2.0.  For  masses  1  and  2  the  effective  b  is 


V  =  51  bt  =  1632  (143) 

bt'  =  1.02  h  =  32.64  (144) 

Equations  (143)  and  (144)  are  similar  to  equations  (130)  and  (131)  for  k.  With  u 
given  by  equation  (133),  Fi  -  1,  and  Ft  -  +1,  equation  (142)  gives  the  amplitudes 
shown  in  equations  (138)  and  (139)  namely 


_ _ 

1632  x  638.7487769  i 


9.5928966 x 10-7  i 
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_ 1 _ 

32.64  x  638.7487769  « 


-  4.7964483  x  IQ"5  » 


Another  test  has  been  run  with  fci  =  0,  b*  =  32,  F\  =  12.5,  Ft  =  .25, 
A;i  =  500  and  kj  =  10s.  Setting  the  force  equal  to  the  mass  gives  an  equal 

acceleration  to  each  body  such  as  that  caused  by  a  gravity  field  variation.  The 

presence  of  the  spring  fci  causes  the  motions  of  mi  and  mj  to  be  slightly  different. 
The  amplitudes  computed  by  program  TWOMASS  at  the  second  resonance  u  = 
638.7493909  rad  /sec  are 

Ai  =  -  2.45121128  x  10-6  -  2.3055  xlO'11  i  (145) 

Aj  =  2.45121128  xl0~6  +  1.15265x10-®  i  (146) 


The  magnitudes  of  A\  and  At  are  essentially  the  real  component.  The  phases 
are 


01  =  -  179.9994610  deg  (147) 

02  =  +  179.9730572  deg  (148) 


Numerical  integration  shows  the  phase  of  zj  to  be  delayed  by  .027  deg  as  expected. 
The  first  maximum  of  z\  is  at  .00491833  sec  and  the  first  maximum  of  ti  is  at 
.00491908  sec. 

A  number  of  other  test  runs  have  been  done  to  check  the  agreement 
between  the  theoretical  expressions  and  the  numerical  integration.  In  the  next 
run  &2  =  016  which  is  t  Q  of  10*  for  mass  2.  The  other  parameters  are  the 
same,  namely  =  0,  =  500,  =  10s.  =  12.5,  and  /:  =  .25.  At  the 

second  resonance,  un  -  638.7493909  rad/sec,  the  theoretically  computed  ampli- 


tudes  are 
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A\  =  -  .245121128  x  10"8  -  .4611069  x  10*7  »  (149) 

At  =  -  .245121128 x 10“8  +  .2305308  x 10" 8  «'  (150) 

The  phase  angles  are 

»i  =  -  178.9223137  deg 
=  136.7569631  deg 

The  simulation  uses  the  initial  conditions 

xi  -  xt  *  -  .245121128  x  10-8  (153) 

*1  =  2.94531779  x  10-8  (154) 

it  =  -  1.4725145x  10-*  (155) 

The  computed  magnitudes  of  Ai  and  At  are  -24516449xl0'i  and  .3364949x10'*. 
The  values  in  the  numerical  simulation  are  .245164x10*  and  .3364949x10"*.  The 
case  was  rerun  with  bt  -  0  and  the  same  initial  conditions.  The  amplitude  of  A\ 
remained  at  .245164  x  10“*  and  the  amplitude  of  At  showed  values  of 

-  .336501  x  10*.  +  .336526  x  10“8,  -  .336551  x  10"*,  +  .336577  x  10~*, 

-  .336602  x  10-*,  and  +  .336627  x  10~8  at  times  .00118165,  .00610012,  .01101859, 
.01593707,  .2085555,  and  .02577403  seconds.  Without  damping,  the  amplitude 
builds  up  at  resonance.  The  next  test  run  was  done  with  bi  =  16  and  bt  =  0 
with  the  other  parameters  kept  the  same  For  this  case  the  computed  amplitudes 


(151) 

(152) 


are 
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Ai  =  -  .2237  x  10-“  -  .11991086  x  10-c  i  (156) 

At  =  -  .1249877  x  10"3  +  .5994955  x  10“*  «  (157) 

at  the  second  resonance,  u>  =  638.74939091  radians/ sec.  The  numerical  simula¬ 
tion  agrees  with  the  theoretical  computation.  In  the  last  simulation,  equal  values 
of  m,  k,  and  b  are  used  for  both  masses  namely, 


mi  r  mj  =  .25 

(158) 

fci  =  bt  =  32 

(159) 

kx  =  kt  =  10* 

(160) 

The  value  of  b  is  approximately  10%  of  critical  damping  for  a  one  mass  system. 
The  resonant  frequencies  are  387 .3580658  and  1032.63630037  rad /sec.  The 
amplitudes  are 

Ax  -  -  .83367611  x  10~6  -  .83992923  x  10~6  i 

Ai  =  -  .1250514  x  10- s  +  .4215206  x  10-6  f  (161) 

at  the  second  resonance.  The  magnitudes  of  Ax  and  At  are  .1183459x10**  and 
.13196459x10**.  The  phases  an  -134.7859257  and  161.3721777  deg.  Figure  20(a) 

shows  i]  as  a  function  of  time  and  Figure  20(b)  shows  rj.  The  first  peak  of  z\ 

* 

is  .1183425x10"*  at  .00227812  seconds.  Since  the  period  is  .006084606  sec,  the 
phase  delay  is  -134.786  deg.  The  first  positive  peak  of  is  .1319643x10"*  at 
.00335715  sec  which  is  a  phase  lag  of  -198.628  deg  or  +161.372  deg.  The  phase  is 
measured  relative  to  the  driving  function  which  is  .25  cos  u >t. 
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The  test  runs  with  the  numerical  integration  program  confirm  the  analytical 
expressions  derived  in  Section  5.0.  The  analytic  expressions  give  only  the  final 
equilibrium  state  whereas  the  numerical  integration  allows  the  transient  behavior 
to  be  studied  also.  In  the  case  of  a  high  Q  system  it  would  be  time  consuming  to 
determine  the  equilibrium  response  by  numerical  integration.  Plotting  resonance 
curves  is  much  more  efficient  with  the  analytic  expressions. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  experiment  configuration  is  now  nearly  frozen  and  the  measurement 
schedule  approaches  crystallization.  The  chances  of  substantial  delay  are 
reducing  every  week  that  goes  by.  We  have  to  fear,  still,  administrative 
delays  in  subcontract  award,  such  as  contract  for  added  scope  to  SAO,  and 
the  contract  for  the  extra  rotating  shutter  to  University  of  Maryland.  On 
the  technical  front,  a  possible  delay  might  have  to  be  faced  in  the  procure¬ 
ment  of  the  custom-made  4°K  cryostat  for  the  cryogenic  force  sensor.  We 
still  have  to  receive  a  quotation  for  it.  Everything  else  is  pretty  much 
under  control . 

Of  great  help  would  be  receiving  from  DARPA  as  soon  as  feasible  the 
full  data  on  the  tritium  source,  so  that  the  high-speed  rotating  table  could 
be  completed,  design-wise,  and  SAO  could  proceed  to  the  issuance  of  purchase 
orders  for  components  and  subsystems.  This  is  particularly  important  becuase 
the  SAO  design  team  is  scheduled  to  go  to  Australia  on  March  15,  and  if  we 
do  rot  have  the  rotating  table  fully  defined  by  then,  there  will  be  a  delay 
of  aDout  six  weeks  in  the  issuance  by  SAO  of  these  purchase  orders  (the  SAO 
team  is  scheduled  to  launch  a  balloon  from  down  under,  with  a  payload  of 
the.r  own  design  and  mechanization). 

Last  but  not  least,  we  recommend  that  DARPA  seriously  consider  funding 
expeditiously  the  added  scope  work. 
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EXECUTIVE  SUMMARY 


During  tne  report  period,  the  project  status  has  evolved  substantially, 
areas  of  uncertainty  have  been  cleared,  and  several  milestones  have  been 
met.  Specific  accomplishments  are: 

(a)  The  torsion  balance  has  been  constructed  and  testeu  at 
University  of  Maryland.  Prof.  Weber  is  awaiting  shipping  instructions 
based  on  the  experiment  location; 

(b)  Same  as  above,  for  the  tuning  fork; 

(c)  The  construction  of  the  cryogenic  force  sensor  has  been  completed 
in  its  room-temperature  version,  and  the  instrument  is  undergoing  laboratory 
testing  at  IFSI-CNR,  Frascati,  Italy.  It  will  be  ready  for  shipment  to 
Raytheon  in  late  August  1988.  In  Portsmouth,  the  instrument  will  be  housed 
in  a  4 6K  cryostat,  system  tested  and  shipped  to  the  test  site; 

(d)  The  1  RPM  rotating  table  for  use  in  testing  the  torsion  balance  is 
now  being  assembled  at  Raytheon,  Portsmouth,  Rl.  It  will  be  ready  for 
shipment  to  the  test  site  in  early  July  1988; 

(e)  A  low-speed  data  acquisiton  system  (suitable  only  for  the  torsion 
Dalance)  has  been  identified  as  Raytheon  available  capital  equipment  and  could 
oe  provisionally  used  at  the  test  site,  starting  in  early  July; 

(f)  A  high-speed  data  acquisition  system  (suitable  for  use  with  all 
three  sensors)  is  now  being  procured  by  Raytheon  and  will  be  ready  for  shipment 
to  the  test  site  in  late  July  1988; 
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(g)  The  problems  (of  both  security  and  safety  nature)  concerning  the 
aelivery  t*'  Raytheon  of  the  100  Kilocurie  source,  would  be  eliminated,  should 
UARPA  choose  Los  Alamos  National  Laboratory  (LAND  as  the  experimental  test 
site.  Sources  of  adequate  strength,  contained  in  enclosures  of  acceptable 
size  and  shape,  are  already  available  at  LANAL-WX-5. 

(h)  With  the  expectation  that  LANL  will  be  chosen  as  the  test  site, 
Raytheon  has  already  contacted  LANL-WX-5,  to  explore  the  possibility  of 
conducting  the  tests  in  a  LANL  lab  area  wh  ich  would  be  assigned  to  the 
project.  Also  LANL  was  asked  when  they  could  provide  Raytheon,  on  a  loan 
basis,  general-purpose  laboratory  instrumentation  and  supplies  needed  during 
testing.  This  includes  the  following  items: 

1.  Multichannel  digital  tape  recorder,  with  A/D  converter  (and  a 
supply  of  blank  tape  reels); 

2.  Digital  clock,  IRIG-code  time  generator; 

3.  Multichannel  chart  recorder  (and  a  supply  of  chart  rolls); 

4.  Liquid  helium  conta iner  and  ref  ills; 

5.  Personal  Computer  PC-AT,  for  use  with  real-time  data  acquisition 
system. 

(i)  Concerning  the  100  Hz  modulator  of  the  neutrino  beam,  for  the 
cryogenic  force  sensor,  Raytheon  could  proceed  by  awarding  to  University  of 
Maryland  a  subcontract  for  the  construction  of  a  large-size  rotating  chopper, 
complete  with  sapphire  crystals.  However,  DAPRA  authorization  for  this 
procurement  has  not  yet  been  received  by  Raytheon  at  the  time  of  this  report. 


An  alternative  way  of  proceeding  would  be  to  have  LANL-WX-5  construct 
(unoer  a  MIPR  from  DARPA)  the  2000  RPM  rotating  table.  A  preliminary 
design  has  been  defined  at  SAO,  as  part  of  Raytheon's  subcontract  to  the 
Observatory. 

(j)  A  computer  simulation  of  the  response  of  all  sensors  to  the 
rotating  source/dummies ,  has  been  initiated  at  SAO,  taking  into  account  the 
dynamic  properties  of  the  torsion  balance,  the  tuning  fork  and  the  cryogenic 
force  sensor.  This  simulation  Is  expected  to  provide  vital  data  concerning 
manufacturing  tolerances  (total  mass  and  position  of  the  center  of  mass) 
required  for  constructing  sources  and  dummies,  to  minimize  the  spurious 
signal  due  to  newtonian  gravity  forces. 

Raytheon  is  planning  a  visit  to  LANL  in  early  June  1988,  to  illustrate  to 
wX-5  the  status  of  the  project,  to  survey  the  laboratory  area  that  could  be 
assigned  to  the  tests  (using  a  100  Kilocurie  tritium  source  in  a  3.5"  sphere, 
or  in  a  similar  container)  and  to  outline  the  measurement  plan  that  Raytheon 
has  formulated.  In  order  to  prepare  for  this  important  visit,  Raytheon  has 
worKed  out  an  updated  presentation,  that  constitutes  the  majority  of  this 
report. 
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RAYTIEON  NEUTRINO  DETECTORS  PROGRW  ORGANIZATION 

WtflER  VERIFICATION  EFFORT 


Hr.  G.  Nystrom 


MEASUREMENTS  TO  BE  PERFORMED  TO  SATISFY  REQUIREMENTS 
OF  CONTRACT  F49620-87-C-0050 
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Large-size  rotating  shutter  (  to  be  built  by  U.  of  Maryland) 

Tripod  with  vibration-abatement  provisions  (  under  construction  at  SAO) 


measurements  to  be  performed  to  satisfy  requirements 
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RAYTHEON  NEUTRINO  DETECTORS  TROCRAM  -  TASK  DESCRIPTION  AND  SCHEDULE  (Continued) 


Data  Collection,  reduction,  process¬ 
ing  and  real-time  analysis 
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POSSIBLE  SOURCES  OF  ?tASUREMENT  ERRORS  THAT  HAVE  BEEN  IDENTIFIED 
AM)  THAT  REQUIRE  SPECIAL  ATTENTION 
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t£B£R  TORSION  BALANCE  -  ELECTRONICS 
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Schematic  diagram  of  overall  instrument 


WO  TUNING  FCRK  AtC  ROTATING  CHOPPER 
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Minimiw  detectable  force  10"'  DYNES 
Integration  time  /v  10?  SEC 
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THE  CRYOGENIC  FORCE  SENSOR  (IN  A  M°K  CRYOSTAT) 
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Minimum  detectable  gradient  (baseline  =  20  cm):  2.83  x  10*3  Eotvos  units 


LSTIMATE  OF  VIBRATION  NOISE 


AN  APPROACH  TO  ABATE  MECHANICAL  VIBRATIONS 
AM)  SEISMIC  OSCILLATIONS 
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By  fiber  optics  visual  observations  from  the 
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BASIC  PLAN  OF  MLASUREMLNTS  AT  LANL 
(A)  Measurements  to  be  carri ed-out  without  the  use  of  a  tritium  source 
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CONCLUSIONS  AND  RECOftENDATIONS 
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EXECUTIVE  SUMMARY 


In  this  quarter,  a  substantial  milestone  was  achieved:  Raytheon  established  a  laboratory 
base  at  LANL,  Los  Alamos,  N.M.  and  received  from  LANL  three  dummies,  for  installation  on 
the  1  RPM  rotating  table.  Specific  accomplishments  were  as  follows. 

1.  Raytheon  1  RPM  rotating  table  and  U.  of  Maryland  room- temperature  Torsion 
Balance  were  installed  in  the  LANL  lab  area  assigned  to  the  project.  The  area  is  in 
Building  86,  Technical  Area  TA-33. 

2.  The  two  pieces  of  equipment  that  LANL  is  constructing  for  our  project  under  a 
direct  DARPA  contract,  (the  support-stand  for  the  1  RPM  rotating  table,  and  the  40" 
table-top,  where  neutrino  sources  and  dummies  will  be  mounted),  are  scheduled  for  delivery 
to  TA-33,  Bldg.  86  in  early  September  88. 

3.  The  low-speed,  PC-controlled,  data  acquisition  system  (based  on  a  Lab/Tech 
Notebook  approach)  is  expected  to  be  ready  for  shipment  by  Raytheon  to  LANL  no  later  than 
September  9,  1988. 

4.  If  the  LANL  neutrino  sources  are  available.  Prof.  Weber's  torsion  balance 
experiment  (with  real-time  presentation  of  the  experiment  results)  could  start  around 
September  12,  1988. 

5.  Raytheon's  cryogenic  force  sensor,  based  on  the  room- temperature  gravity 
gradiometer  developed  at  SAO/IFSI-CNR  by  Co-Investigator  Dr.  Fuligni  is  expected  to  be 
ready  for  4°K  laboratory  tests  in  Portsmouth,  RI  in  the  last  week  of  September  1988. 
According  to  the  DARPA  contract,  these  tests  will  not  require  neutrino  sources  and  will 
be  based  on  the  measurement  of  instrument  sensitivity  to  electrostatically  generated 
"force  signals".  High-speed  data  collection  (the  cryogenic  force  sensor  resonates  at 
100  Hz,  while  Weber's  torsion  balance  has  a  resonance  period  of  16  seconds)  will  be 
performed  with  the  Diabase  approach.  The  software  for  Diabase  does  not  function  yet 
and  will  require,  most  probably,  the  assistance  of  the  software  developers,  Morra  Inc. 

The  above  instrumentation  will  be  shipped  to  Los  Alamos  and  tested  there  with 
neutrino  sources  when  the  2000  RPM  rotating  table  is  available.  DARPA  has  not  yet 
funded  the  related  development  effort.  A  possible  date  of  availability  of  this  table 
is  January  1989,  at  the  earliest. 


6.  The  SQUID  magnetometer  (built  under  a  1988  Raytheon  in-house  IDP  project)  was 
successfully  assembled  and  tested  in  Portsmouth.  Based  or.  DARPA's  contract,  we  are 
scheduling  tests  at  Raytheon  on  basic  feasibility  issues  concerning  the  potential  use 
of  the  SQUID  in  detecting  neutrinos,for  the  second  half  of  September  1988.  As 
contractually  specified,  these  tests  will  not  require  neutrino  sources  and  will  consist 
of  the  following  tasks:  (a)  measurement  of  minimum-achievable  SQUID  intrinsic  noise, 
when  the  SQUID  is  loaded  by  a  large-mass,  high-permeability  interaction  target;  (b) 
measurement  of  maximum-achievable  collecting  area,  using  superconducting  transformer; 

(c)  measurement  of  maximum- achievable  relative  permeability  for  interaction  target. 

As  can  be  seen  from  items  1  through  6  above,  Raytheon  effort  concentrated  in  getting 
ready  for  the  verification  of  the  results  that  Prof.  Weber  claimed  he  has  obtained  with 
his  room- temperature  torsion  balance.  This  was  the  focus  of  project  activity  in  the 
past  quarter.  Progress  was,  however,  also  made  in  the  area  of  the  cryogenic  force  sensor. 
This  will  become  the  next  focus  of  activity,  once  the  tests  on  the  torsion  balance  have 
been  carried  out.  Also,  we  have  moved  further  in  another  area:  the  development  of  the 
magnetic  interaction  sensor.  For  this  last  instrument,  DARPA  has  only  authorized  laboratory 
tests  on  basic  (and  controversial) ,  feasibility  issues  pertaining  to  magnetic  material 
properties  and  high-sensitivity  SQUID  development,  with  the  proviso  that  Raytheon  construct 
with  its  own  IDP  funds,  as  an  in-house  effort,  the  SQUID  itself  (which  Raytheon  did). 

By  the  time  the  next  quarterly  report  is  submitted  (15  November  1988),  we  expect  that 
final  conclusions  will  be  available  concerning  Weber’s  torsion  balance,  thus  completing 
the  most  immediate  and  urgent  assignment  that  Raytheon  received  from  DARPA. 


iv 


STATUS  OF  READINESS  FOR  INSTRUMENTATION  TO  BE  USED  AT  LANL,  LOS  ALAMOS,  N.M. 


1.  Room- temperature  Torsion  Balance 

The  room- temperature  torsion  balance,  that  Prof.  Weber  has  constructed  at  U.  of 
Maryland  under  the  terms  of  a  Raytheon  subcontract,  was  installed  in  late  July  1988 
at  LANL.  Figures  1,  2  and  3  show  the  mounting  of  the  instrument.  The  figures  depict 
a  Vac-Ion  cryogenic  pump  which  allows  torsion  balance  measurements  in  vacuum.  At 
this  point,  the  torsion  balance  has  been  operated  at  LANL  without  vacuum  by  Prof. 

Weber  and  his  assistants.  With  the  system  addition  of  an  electronic  damper,  the 
balance  exhibited  a  damping  time  of  about  39  seconds,  and  a  resonance  period  of  about 
16  seconds.  The  torsion  balance  was  tested  by  Prof.  Weber,  after  installation  at 
TA-33,  LANL  and  it  was  found  to  meet  pre-shipment  performance.  Weber  made  a  crude 
evaluation  of  sensitivity  by  moving  two  lead  bricks  near  the  balance's  disk  and 
seeing  that  the  response  was  identical  to  what  he  obtained  in  Maryland  with  a  similar 
arrangement.  There  was  also  no  change  in  the  tendency  of  the  balance  to  drift  away 
from  a  mean  position.  Because  of  this  drift,  Raytheon  designed  its  data  collection 
system  with  the  ability  to  search  for  a  signal,  smaller  than  1  millivolt,  on  an 
overall  output  that  might  be  a  large  as  ±  10  volts.  As  requested  by  Prof.  Weber, 
Raytheon  Analog-to-Digital  converter  will  have  a  resolution  of  16  bits  (instead  of 
the  12  bits  used  by  Prof.  Weber  in  his  tests  at  U.  of  Maryland).  As  a  consequence, 
the  overall  instrumentation  system  that  Raytheon  is  assembling  at  LANL,  TA-33,  will 
have  superior  performance  than  that  achieved  in  Prof.  Weber's  torsion  balance  tests 
at  U.  of  Maryland. 

2.  1  RPM  Rotating  Table 

Figure  4  and  5  show  the  rotating  table  that  Raytheon  has  installed  at  TA-33.  There 
are  still  two  pieces  missing  from  the  table:  the  support-stand  and  the  table- top. 

LANL  will  construct  these  under  the  terms  of  a  direct  DARPA  contract.  The  performance 
characteristics  of  the  table  are  given  in  Table  1.  The  addition  of  the  two  missing 
subsystem  will  not  affect  this  performance.  The  table  is  controlled  from  a  H-P 
computer  master  control,  visible  in  Figure  4. 
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Figure  1- 

Front  view  of 
Torsion  Balance 


Figure  2  - 
Another  front 
view  of  the 


torsion  balance 


■jr%, 


Figure  5  - 
The  1  RPM 
rotating  table 


TABLE  1 

PERFORMANCE  SPECIFICATIONS  OF  ROTATING  TABLE 
MODEL  SA  5115B-3 

1.  Table  out  of  Roundness  . 0.015"  3.81  10  * 

2.  Caliper  out  of  Roundness  . 0.010"  2.54  10  * 

3  - -  Ax1s  . 0.005"  . 1.27  10”1 

4.  Flatness  of  Top . 0.0005"  1.27  10”* 

5.  Rocking  . 0.003"  7.62  10”* 

(*)  This  could  increase  with  addition  of  table- top 


mm  (w/dial  meter) 
mm 


mm 


mm 


mm  (*)  Outside 
Diameter 


3.  SQUID  Magnetometer  at  Raytheon,  Portsmouth,  RI  for  Tests  on  Magnetic  Interaction  Targets 

During  this  past  quarter,  Raytheon  has  assembled  and  tested  a  SQUID  magnetometer,  as  a 
part  of  the  company's  1988  IDP  Program.  Next  quarter,  we  plan  to  use  this  SQUID  in  critical 
feasibility  tests  of  direct  interest  to  DARPA's  neutrino  detectors  program,  funded  by  the 
DARPA/AF0SR  contract.  The  magnetometer  consists  of  basically  four  units:  the  SQUID  proper, 
the  cryogenic  descendent,  the  preamplifier  head  with  filters,  and  finally  the  control  unit. 
The  preamplifier  head  is  connected  to  the  cryogenic  descendent.  At  the  other  (lower)  end  of 
the  descendent,  inside  a  superconducting  magnetic  shield,  is  the  SQUID  proper.  The  control 
unit  is  energized  by  rechargeable  batteries  with  about  a  10  hour  lifetime.  We  have  conducted 
lab  tests  to  measure  the  intrinsic  noise  of  the  system  and  found  it  to  be  about  1.8  10"® 
d/H z*5,  where  ff  s  2.07  10’1®  Weber,  for  frequencies  above  10  Hz.  Commercial  SQUIDs  are 
characterized  by  an  intrinsic  noise  of  1.5  10  /  /Hz  ,  in  comparable  conditions. 

As  far  as  the  energy  sensitivity  is  concerned,  our  SQUID  is  characterized  by  a  value  of 
8.5  x  10”31  J/Hz,  against  a  value  of  3  x  10’3^  J/Hz  typical  of  commercial  SQUIDs.  There 
is  therefore  an  improvement  of  about  18  db  in  intrinsic  noise  and  of  about  6  db  in  energy 
sensitivity,  when  compared  to  commercial  units.  We  are  now  in  a  position  to  perform  the 
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three  sets  of  measurements  to  address  the  feasibility  of  the  magnetic  interaction  sensor: 

(a)  measurement  of  minimum  achievable  intrinsic  noise,  when  loading  the  SQUIO  with  a 
heavy,  high-permeability  magnetic  interaction  target;  (b)  measurement  of  maximum  achievable 
collecting  area,  using  a  superconducting  transformer;  (c)  measurement  of  maximum-achievable 
relative  permeability  of  the  target  materials.  We  plan  to  perform  these  measurements  in 
the  forthcoming  quarter  of  contract  performance. 

4.  The  Cryogenic  Force  Sensor 

The  fundamental  subsystem  of  this  instrument  is  the  room- temperature  gravity  gradiometer 
nearing  completion  at  IFSI-CNR,  Frascati,  Italy.  This  subsystem  is  expected  to  reach 
Raytheon  Portsmouth,  RI  in  the  third  week  of  September  1988.  The  4°K  cryostat  will  reach 
Raytheon  in  the  same  time  period,  so  that  by  the  end  of  FY  88,  the  cryogenic  force  sensor  will  be 
undergoing  system  tests  by  the  use  of  electrostatically-generated  "force-signals".  Transfer 
of  this  instrument  to  LANL,  TA-33  and  testing  it  with  neutrino  sources  must  wait  for  DARPA's 
decision  concerning  the  construction  of  the  2000  RPM  rotating  table.  As  contractually 
allowed,  Raytheon  will  keep  this  instrument  here  in  Portsmouth,  RI,  until  the  high  speed 
table  is  available. 


-6- 


CONCLUSIONS  AND  RECOMMENDATIONS 


No  major  obstacles  are  expected  to  arise  in  Raytheon's  feasibility  testing  of  detecting 
low-energy  neutrinos  with  Prof.  Weber's  room- temperature  torsion  balance.  Concerning 
possible  tests  with  the  cryogenic  force  sensor,  developed  by  SAO/IFSI-CNR  under  a  Raytheon 
subcontract,  everything  has  proceeded  according  to  plans,  apart  from  a  few  weeks  delay. 
However,  final  neutrino  tests  cannot  be  scheduled  and  performed  until  the  2000  RPM  rotating 
table  is  available  at  LANL.  Since  April  8fi,  DARPA  has  had  necessary  information  to  make 
a  decision  to  fund  this  effort.  We  recommend  that  this  decision  be  made  without  further 
delay.  Finally,  as  far  as  the  magnetic  sensor  is  concerned,  we  have  to  perform  the  basic 
feasibility  tests  on  the  SQUID  system  and  magnetic  materials,  before  we  can  hope  to  receive 
DARPA  authorization  to  proceed  on  the  procurement  of  the  large-size  interaction  target 
and  the  large  4°K  cryostat  required  by  this  instrument.  Any  decision  on  this  must  wait 
for  the  results  of  the  on-going  tests. 
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-  miwj*  =  lO4^!  (91) 

Therefore  from  equation  (88)  the  motion  of  mass  mi  is  still  given  by  equation  (89) 
to  an  accuracy  of  about  1  part  in  104.  This  provides  a  qualitative  explanation  of 
the  factor  of  104  seen  between  cases  1  and  S. 

In  case  4,  the  high  Q  of  system  2  is  somewhat  degraded  by  the  low  Q  of 
system  1.  Since  there  is  a  factor  of  10*  between  the  resonant  frequencies,  there  is 
not  a  large  damping  effect  on  system  2  because  of  the  low  Q  in  system  1. 

The  calculations  in  cases  1  to  4  were  all  done  at  the  u>2  resonance  computed 
by  setting  equation  (56)  equal  to  sero.  For  cases  1,  2,  and  3,  we  have 

=  638.7493909  rad/sec  (92) 

For  case  4  we  have 

u*  =  638.7494069  (93) 

In  the  presence  of  damping  the  maximum  amplitude  may  not  coincide  exactly 
with  u>:  calculated  from  equation  (56).  In  order  to  determine  the  actual  positions 
of  the  maximum  amplitude,  resonance  curves  have  been  computed  using  the  full 
solutions  given  by  equations  (76)  and  (77).  For  case  1,  the  actual  position  of  the 
maximum  of  A\  is  <*>:  =  638.7493893  rad/sec  and  for  case  4,  u>:  =  638.7493874 
rad/sec.  For  case  1  the  frequency  shift  is  >.0000016  rad/sec.  For  case  4  it  is  : 
.0000195  rad/sec.  Plots  of  Ar  show  it  to  be  sero  at  exactly  the  values  given  in 
equations  (92)  and  (93).  In  both  cases  the  frequency  shift  is  negligible  compared 


